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This article reviews the quantum entanglement in Valence-Bond-Solid (VBS) states 
defined on a lattice or a graph. The subject is presented in a self-contained and peda- 
gogical way. The VBS state was first introduced in the celebrated paper by I. Affleck, 
T. Kennedy, E. H. Lieb and H. Tasaki (abbreviation AKLT is widely used). It became 
essential in condensed matter physics and quantum information (measurement-based 
quantum computation). Many publications have been devoted to the subject. Recently 
entanglement was studied in the VBS state. In this review we start with the definition of 
a general AKLT spin chain and the construction of VBS ground state. In order to study 
entanglement, a block subsystem is introduced and described by the density matrix. Den- 
sity matrices of 1-dimensional models are diagonalized and the entanglement entropies 
(the von Neumann entropy and Renyi entropy) are calculated. In the large block limit, 
the entropies also approach finite limits. Study of the spectrum of the density matrix led 
to the discovery that the density matrix is proportional to a projector. 

Keywords: AKLT; VBS; Quantum Spins; Entanglement; Entropy; Density Matrix. 



1. Introduction 

The fields of statistical mechanics, condensed matter physics and quantum infor- 
mation theory share a common interest in the study of interacting quantum many 
body systems. The concept of entanglement in quantum mechanics has significant 
importance in all these areas, it was introduced in^^. Roughly speaking, entangle- 
ment is a phenomenon of quantum mechanical nature in which quantum states 
of physical systems are linked together so that one system can not be adequately 
described without full mention of its counterpart, even when the individual sys- 
tems may be spatially separated. Entanglement becomes particularly interesting in 

1 
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a many body interacting system where a subsystem strongly correlates with its en- 
vironment (other parts of the system). The correlations may reject the principle of 
local realism, which states that information about the state of a system can only 
be mediated by interactions in its immediate surroundings (neighbors). The charac- 
teristic length of entanglement may be diverging while the usual correlation length 
remains finite Quantum entanglement is a fundamental measure of 'quantum- 
ness' of a system: how much quantum effects we can observe and use to control one 
quantum system by another. It is the primary resource in er nergin g technologies 
of quantum computation and quantum information processing l ^ l ^^ l. Entanglement 
properties play an important role in condensed matter physics, such as phase transi- 
tions "^211 and macroscopic properties of solids There is an excellent review 1221 
on entanglement detection, which covers many topics (such as multi-partite entan- 
glement, concurrence. Bell inequalities, entanglement witnesses and experiments) 
beyond the scope of our present discussion. 

Much of current research seeks to elucidate quantum entanglement in a 
variety of interacting systems. Extensive research has been undertaken to in- 
vestigate quantum entanglement in strongly correlated states such as spin 
chains, corre lated electrons, interacting bosons as well as othe r models. (See 
I3I4I9I11I31I151 23 24 38 41 42 44 46 47 50 55 56 57 62 65 66 67 74 7W9l ^j. reviews and 

references.) A general approach studies the density matrix of a certain subsystem 
of a strongly entangled state. Th e co ncept of the reduced density matrix was first 
introduced by P. A. M. Dirac in 1131. Following C. H. Bennett, H. J. Bernstein, S. 
Popescu and B. Schumacher 1^, the entropy (derived from the spectrum of the den- 
sity matrix) of a subsystem serves as the measure of entanglement for a pure state. 
The von Neumann entropy and its generalization (the Renyi entropy) are typical 
quantifications of entanglement. These characteristic functions may depend on the 
physical parameters (size, coupling constants, external fields, etc.) in various differ- 
ent ways. An area law for t he von Ne umann entropy in harmonic lattice systems 
has been extensively studied | 12 | 13 | 40 | j^^ie area law states that the entropy scales 
proportional to the size (area) of the boundary of the subsystem (this statement was 
strictly proved for gapped models only) . The entropy of the whole system vanishes 
if the system is in a pure state (usually the unique ground state), but it can be pos- 
itive for a subsystem (this means entanglement). This is in contrast with a classical 
system in which if the entropy is equal to zero for the whole system, then it also 
vanishes for any subsystem. In quantum mechanics, when the whole system is in a 
pure state, a subsystem can be in a mixed state, e.g. The Bell pair (or the Einstein- 
Podolsky- Rosen state) (| — | |t))/\/2 is a pure state. A single spin subsystem 
is in a completely mixed state. The density matrix of one spin in an EPR state is 
proportional to identity matrix: p = (| t)(t I + I l)/2, which has entropy In 2. In 
a more general case, we shall consider a pure system consisting of two subsystems 
B (Block) and E (Environment) with two different Hilbert spaces Hb and of 
dimensions diniB and dimE, respectively (assuming that dims < diniE). So the 
entropy of the whole system B U E is zero. If the entropy of a subsystem (B) is 
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positive then the wave function of the whole system does not factorize: 

\BUE)^\B)®\E). (1) 

Let and {|efe)} be any fixed orthonormal bases for subsystems B and E, 

respectively. Then the wave function of the whole system can be written as 

|i3Ui?)-^A,fc|6,)0|e,), (2) 

jk 

for some matrix A of complex elements Ajk- Then by the singular value decompo- 
sition (see §2.1.10 oi^, A = UVV, where V is a diagonal matrix with non- 
negative elements, and U and V are unitary matrices {A is a rectangular ma- 
trix if diiriB 7^ dimE, in that case U and V have different dimensions). Define 
\Bi) = J2j Uji\bj), \E^) = J2k Vik\ek), and y/Xl = Va. With such a choice of bases, 
we find that the wave function can be put in the following form: 

D 

|i?Ui?) =^V^|B,)®|£;,) (3) 

where as well as {|£'j)}, are orthonormal bases. The coefficients < ^yXj < 1 

satisfy 

EA,=1. (4) 

Once again, if the entropy of the subsystem S[B] > then D > 1 {D < diniB)- In 
this case we can control one subsystem by the other. Indeed if we measure subsystem 
B (in \Bj) basis), then the wave function of the whole system will change (collapse) 
into 

|S.o) ® \E,,) (5) 

with probability for some jo- So that a measurement of B changes the state of 
E and this is quantum control. This is an important resource for quantum device 
building (and maybe even quantum computation). Much insight in understanding 
entanglement of quantum systems has been obtained by studying exactly solvable 
models in statistical mechanics, in which it is possible to solve the subsystem density 
matrix and calculate the entropy exactly. This review is devoted to one particular 
set of models - the AKLT models. 

In 1987, 1. Affieck, T. Kennedy, E. H. Lie b and H. Tasaki proposed a spin inter- 
acting model known as the AKLT model 1^^. The model consists of spins on lattice 
sites and the Hamiltonian describes interactions between nearest neighbors. The 
Hamiltonian density is a linear combination of projectors. Each projector is written 
as a polyno mial of the inner product of a pair of interacting spin vectors. The authors 
(AKLT) of ^1^2] found the exact ground state, which has an exponentially-decaying 
correlation function and a finite energy gap. In their early works, authors discussed 
the 1-dimensional AKLT lattices with open and periodic boundary conditions and 
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2-diniensional models such as the hexagonal lattice. D. A. Arovas, A. Auerbach and 
F. D. M. Haldane studied the model in the Schwinger boson representation (see 
S I2.4.2I below) and calculated the correlation functions using the coherent state ba- 
sis (see S I5.1.2I below). In their work 13, it was shown that the quantum spins in the 
AKLT model is equivalent to clas sical spins. TM^^ has been attracting enor- 

mous research interests since t hen I _ l _ l | ol | o/ | c)4 | /l | -^^ defined and solved 

in higher dimensional lattices HHnEIIZI] jg generalizable to the inhomogeneous 
(non-translational invariant) ca se (spins at different lattice sites may take different 
values) and an arbitrary g raphi I 1 ^ The AKLT model has also been general- 
ized to the SU(n) version ~j- .ffX.^?^ certain conditions (see §[23]), the 
ground state has proven to be unique HEEIEll, it is known as the Valence-Bond- 
Solid (VBS) state. The VBS state is interesting to different research fields. The 
Schwinger boson representation of the VBS s tate (see (jSHl)) relates to the Laughlin 
ansatz of the fractional quantum Hall effect I 5 | 38 | 43 | ^]^^ Laughlin wave function 
of the fractional quantum Hall effect is the VBS state on the complete graphic. 
The VBS state illustrates ground state properties of anti-ferromagnetic integer-spin 
chains with a Haldane gap^^. In 1-dimension, the VBS state is related to the ma- 
trix product state discovered by A. Kliimper, A. Schadschneider and J. Zittartzand, 
ggg [53E11 -y^rhere a 'q-deformed' VBS-model (a generalized anisotropic version) was 
studied. The geometric entanglement for a 1-dimensional spin-1 VBS state was 
studied in EH. The VBS state can be used as a resource state in measurement-based 
quantum computation invented by R. Raussendorf and H. Briegel^^. VBS st ate al- 
lows universal quantum computation 1^ . Many people developed these ideas l^^l^. 
AKLT Hamiltonian can be implemented in optical lattices 

This review introduces some of the main results on quantum entanglement in 
VBS states defined on a lattice or a graph. We take a pedagogical approach, starting 
with the basics of the AKLT model, the construction of VBS states and measures of 
entanglement (special attention being paid to the uniqueness of the ground state for 
a finite lattice or graph). We shall consider a part (subsystem) of the ground state, 
i.e. a block of spins. It is described completely by the reduced density matrix of the 
block, which we call the density matrix later for short. The density matr ix has been 
studied extensively. It contains information of all correlation functions | 5 | 46 | 48 | 7 6l 
The entanglement properties of the VBS stat es has been studied by means of the 
density matrix as inl l'i' l l8 | 19|25|48 | 74 | 76|7 Z178j ^^^^ ^Yi^^^ throughout the review we 

are considering zero temperature so that the AKLT spin system lies in the VBS 
ground state. 

The review is divided into eight sections including a complete treatment of 1- 
dimensional models with open boundary conditions: 

(1) A brief introduction to the topic. (§[1]) 



*'Many more people work on and contribute to the subject and our list of references is far from 
complete. 
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(2) The construction of the general AKLT Hamiltonian. Introduction of the VBS 
ground state. Definition of different versions of the AKLT model: 1) The basic 
model; 2) The generalized (including the inhomogeneous) model. Proof of the 
uniqueness of the VBS ground state on a finite graph. In this and the next 
sections, we consider the AKLT model on an arbitrary connected finite graph. 
This includes all lattices in any dimension. (§[2]) 

(3) In order to study entanglement, the graph (or lattice) is divided into two sub- 
systems (the block and the environment). Different measures of entanglement, 
namely, the von Neumann Entropy and the Renyi entropy are defined. We de- 
fine anther Hamiltonian called the block Hamiltonian. The block Hamiltonian is 
the AKLT Hamiltonian for the block, but the uniqueness condition is violated. 
The block Hamiltonian is used to describe general properties of the block den- 
sity matrix. We discuss the relation between the spectrum of the density matrix 
and the degenerate ground states of the block Hamiltonian. (§[3]) 

(4) The simplest 1-dimensional AKLT model with spin-1. Calculation and diago- 
nalization of the density matrix (in an algebraic approach) . Calculation of the 
entanglement entropies. (§[4]) 

(5) The 1-dimensional AKLT model with spin-S". Calculation and diagonalization of 
the density matrix using the Schwinger boson representation. Discussion of the 
entanglement entropies. Derivation of the relation between the density matrix 
and correlation functions. (§[5]) 

(6) The 1-dimensional inhomogeneous model (spins at different lattice sites can be 
different). Discussion of entanglement entropies. (§0) 

(7) The 1-dimcnsional SU{n) model (in the adjoint representation). Calculation of 
entanglement entropies. (§[7]) 

(8) A summary including open problems and conjectures. (§[S1) 

2. The General AKLT Model 

In the following we give the most general AKLT Hamiltonian and VBS states. The 
definition applies to both graphs and arbitrary lattices. 

2.1. The Hamiltonian 

The original AKLT Hamiltonian describes a spin interacting system, in which spins 
sitting at lattice sites interact with nearest neighbors. One of the most simple ver- 
sions is an (1-dimensional) open chain of sites with spin-1 at each site, and the 



site j. This Hamiltonian ((5]) looks like the Heisenberg Hamiltonian with an extra 



Hamiltonian is given by ^ 




Sj , S^l , S j ^ at 



(6) 
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quadratic term (the proportionality factor 1/2 and the additive constant 2/3 are 
sometimes neglected which only shifts and scales the energy spectrum as a whole), 
but the physical system behaves quite differently. It was later generalized: 1) The 
spin Sj at each site can take higher values; 2) Different lattice sites can have different 
spins; 3) Different boundary conditions {e.g. an open boundary condition) can be 
applied. An arbitrary boundary condition or distribution of spin values over sites 
may not yield a unique ground state {e.g. Hamiltonian ([6]) has 4- fold degenerate 
ground states). We could find the condition for the uniqueness of the ground state 
(see § 12.51 and . The Hamiltonian can be d efine d on higher dimensional lattices 
{e.g. 2-dimensional square or hexagonal lattice or an arbitrary graph (A graph 
consists of two types of elements, namely vertices and edges. Every edge connects 
two vertices. Spins locate on vertices and two vertices connected by an edge are 
nearest neighbors. An edge is also called a bond physically. In our considerations 
only nearest neighbors are present in the Hamiltonian. More formal and detailed 
explanation of a graph is given in ii l2.3.1l ). These different versions of the AKLT 
Hamiltonian under consideration share two common features: 

(1) The Hamiltonian is a sum of terms with only nearest neighbor interactions, i.e. 

H = Y, H{k, 0, {kl) e {edges}. (7) 

(kl) 

Here the Hamiltonian density H{k, I) describes the interaction between two 
spins at vertices k and / for a connected graph. The construction for a lattice 
is similar. Only spins at nearest neighbor pairs (A;Z)'s called bonds (edges of a 
graph) interact. 

(2) The Hamiltonian density H{k, I) is a sum of several terms. Each term is propor- 
tional to a projector. The proportionality coefficients are all positive numbers. 
i.e. 

i?(fc,0=^Cj(fc,07rj(A;,Z). (8) 
J 

Here 7rj(/c,/)'s are projectors and Cj(fc, Z)'s are positive coefficients. Note that 
the coefficients may depend on the bond {kl) and the projector TTj{k,l). The 
projector 7rj(fc, /) projects Sk and Si on the joint value J of the bond spin 
Jki — Sk + Si. (We also call the joint spin value J the bond spin value) . The 
meaning is this: The spin Sk with spin value Sk at site (or vertex) fc is a (25*^+1)- 
dimensional representation of the SU{2) Lie algebra, while Si is a (25*; + 1)- 
dimensional representation. The direct product of these two representations is 
reducible to a direct sum of irreducible representations with dimensions 2 J + 1 
and J runs from — S*;] to Sk+Si. The Hilbert space 'splits' into these invariant 
subspaces labeled by J which is called the bond spin value of S'fc and Si. (The 
eigenvalues of the Casimir operator - the square of the bond spin {Sk + Si)^ is 
J {J + 1).) The projector 7rj(fc, /) projects on the invariant subspace with bond 
spin J. If we choose an orthonormal basis {| J, to) | m = — J, . . . , J} for the 
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subspace, such that {Sk + Sif\J,m) = J(J + l)\J,m) and (5^ + S'f)|J, m) = 
m|J, m), then the projector could be written as 

J 

^j{k,l)^ |J,m)(J,m|. (9) 

m— — J 

This form ^ is cumbersome in practical use and it is preferred to express the 
projector TTj(k,l) explicitly in terms of spin operators Sk and Si. We shall do 
that in the next section fii l2.2l ). 

Even without an explicit form of the projectors, an immediate consequence of these 
two properties is that the Hamiltonian is positive semi-definite. Furthermore, 
because of this, if we could construct a state {tp) which has no projection on any of 
the specified bond spin- J states appearing in ^ for each bond, i.e. 7rj(fc, = 0, 
V (kl), then it has to be a ground state (with energy equal to zero) regardless of the 
specific values of the coefficients. El The uniqueness condition of the ground state 
will be discussed later in § 12.51 . 

2.2. The Projector 

In order to complete the definition of the general AKLT Hamiltonian ([7]) and Hamil- 
tonian density (jS]), we have to give an explicit expression of the projector TTj{k,l) in 
terms of spin operators Sk and Si. We derive the expression in two steps. The forms 
of 7rj(fc, I) for a specific model such as the expression © or those for 1-dimensional 
models in §[4], §[5] and §[6] can all be obtained through this approach as follows. (An 
explicit construction of the projector was given inl^.) 

(1) Consider the following two sets of operators: the projectors 

{nj{k,l) \J=\Sk-Si\,...,Sk + Si} (10) 

and the powers of the inner product {Sk ■ Si) 

{{Sk-Sir \n^0,...,2S<}, 5< =min{5fc,5z}. (11) 

One set is expressible in terms of the other. They are related by a linear trans- 

''The Hamiltonian is essentially a sum of projectors with positive coefficients. A projector n satisfies 
TT^ = n. So tfiat for an arbitrary state we fiave {ip\n\r(>) = {tp\n^\tp) > because it is an inner 
product of 7r|-i/)) with itself. 

'^Some authors add or omit additive constants in the expression of projectors, e.g. in ^ the 
constant 2/3 can be dropped. This may shift the ground state energy but does not affect the form 
of the ground state because the ground state is constructed to have no projection on the specified 
subspaces for every bond. 
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formation: 



Sk ■ Sir 

-j [{Sk + Sif -Sk{Sk + l)-Si{Si + l)]'' ^'^(^'') 

^ J=\Sk-S,\ 

E (9) [JiJ+^)-Sk{Sk + l)-Si{Si + l)r7rj{k,l) 

(12) 



J=\Sk-Si\ 



for n = 0, . . . , 2S'< . In we have used 

Sk+Si 

E ^Jikj) = l (13) 

J=\Sk-Si\ 

being the identity. This set of 25< + 1 hnear equations can be inverted, 
which express the projector 7rj(fc, I) as a polynomial of the inner product {Sk ■ 
Si). We shall not pursue with the inversion of but to construct the projector 
in the next step. 

(2) The next step is to realize that if an operator V{k,l) satisfies the following 
conditions 

V{k,l)TTj,{k,l)^SjJ,TTj{k,l), VJ' (14) 

then the operator V{k, I) is identified with 7rj(fc, /) because 

Sk+S, 

r{kj)^r{k,i) J2 '^■^'(^'0 

J'=\Sk-Si\ 
Sk+Si \ 
J2 SjJ,\TTjik,l)=7Tj{kJ). (15) 

J'=\Sk-Si\ / 

Therefore we could construct an operator satisfying the condition (jl4p : 

\Sk-Si\<]<Sk+Si ^ J J\J J 

When 'P{k, I) acting on Trji{k, I), we have: 

(a) If J' ^ J, then the numerator of one factor in the product vanishes, so that 
V vanishes, i.e. V{k,l)nj'{k,l) — 0, if J' ^ J; 

(b) If J' ~ J, all factors in the product become equal to 1, so as the expression 
V. i.e. V(k,l)Trj{k,l) = iTj{k,l). 

So that (fT6|) is the projector 7Tj{k,l), i.e. 7Tj{k,l) — V{k,l). Operator 

... rr {Su + sif-3{j + i) 

n j(j + i)-j(^- + i) ^''^ 

\Sk-Si\<j<Sk+Si ^ J J\J I 
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projects the bond spin J^i = Sk + Si on the subspace with fixed total spin value 
J and \Sk-Si\ < J < Sk + Si. Note that we could expand (Sk + Si)^ = 2Sk-Si + 
•S'fc ( 5fc + 1 ) + 5; ( 5; + 1 ) . Therefore pro j actor TT J ( fc , / ) in (flT)) is a polynomial in the 
scalar product (Sk - Si) of degree 25'<, where S'< = niin{S'fe, Si} is the minimum 
of the two spin values of the same bond (kl). For example with Sk — Si — 1, 
we may have a quadratic polynomial as in (jB)): 

7r2(fc, = ^{Sk ■ Sif + ^{Sk ■ Si) + i (18) 

2.3. The Basic AKLT Model 
2.3.1. The Hamiltonian 

Once we have the building blocks for the Hamiltonian from ? I2.1I and ^ \2.2\ . various 
types of the AKLT model can be constructed. Let us start with the definition of the 
basic AKLT model on a connected graph or lattice. (Any lattice is a special graph 
with periodic structure; our notations and definitions refer to the most general). 
A graph consists of two types of elements, namely vertices and edges. Every edge 
connects two vertices. As in Figure ([T]), a vertex is drawn as a (large) circle and an 
edge is drawn as a solid line connecting two vertices. For every pair of vertices in 
the connected graph, there is a walk El from one to the other. Vertices can also be 
called sites and edges sometimes called links or bonds. In a graph, a pair of vertices 
connected by an edge is regarded nearest neighbors, i.e. the terms edge, bond and 
nearest neighbor are equivalent and interchangeable. (For a lattice, vertices become 
sites and bonds become lattice vectors connecting nearest neighboring sites.) In the 
case of a disconnected graph, the Hamiltonian ([7]) is a direct sum with respect to 
connected components and the ground state is a direct product. Therefore we shall 
need only to study a connected graph. Also, assuming that the graph consists of 
more than one vertices to avoid the trivial case where there would be no interaction 
at all. 

Let us introduce mathematical notations. By Si we shall denote the spin op- 
erator located at vertex I with spin value Si. In the basic model we require that 
Si — zi/2, where zi is the number of incident edges connected to vertex Z, also 
known as the valence or coordination number (the number of nearest neighbors of 
the vertex I). The relation between the spin value and coordination number must 
be true for any vertex I, including boundaries. This will guarantee the uniqueness 
of the ground state, see ? I2.5I . For a lattice, this condition would also yield that bulk 
spins (spins not on the boundary) take the same value z/2 because the number of 
nearest neighbors z is a constant. 

In the basic model we define the Hamiltonian density H{k, I) for bond (edge) 



'^A walk is an alternating sequence of vertices and edges, beginning and ending with a vertex, in 
which each vertex is incident to the two edges that precede and follow it in the sequence, and the 
vertices that precede and follow an edge are the endvertices of that edge. 
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{kl) as 

H{k, I) = Cik, I) 7rs,+s, {k, I), H{k, I) > (19) 

with C(fc, I) an arbitrary positive real coefficient (it may depend on the bond (kl)). 
So that the Hamiltonian density for each bond (edge) is proportional to the projector 
on the subspace with the highest possible bond spin value {Sk + Si). The physical 
meaning is that interacting spins do not form the highest possible bond spin (this 
will increase the energy) in the ground state. Then the basic AKLT Hamiltonian 
on an arbitrary connected graph according to ((T]) is 

H^Y. = E ^(^' ^s,+s, {k, I). (20) 

(kl) (kl) 

Here we sum over all bonds (edges) {kl) . Note that for a lattice all the highest bond 
spin values for bulk spins are the same and equal to the coordination number z. 
For example, the basic model defined on a 2-dimensional square lattice must have 
spin-2 in each bulk vertex and t:^ in the Hamiltonian; also, the Hexagonal lattice 
has spin-3/2 in each vertex and tts in the Hamiltonian. 

2.3.2. The VBS State - Pictorial Method 

In this section we consider the Hamiltonian (|20p and construct a ground state which 
is denoted by |VBS). Later we shall see in §[2T5] that it is the unique ground state. 
The Hamiltonian ([20)1 with condition 

Si = \zi (21) 

has a unique ground state | 1 | ^ I 5 | 52 | ]jjiown as the Valence-Bond-Solid (VBS) state. 
It can be constructed in a Pictorial way as follows (see Figure [ij . Each vertex I 
is represented by zi spin-l/2's. We associate each spin-1/2 with an incident edge. 
In such a way each edge has two spin-l/2's at its ends. We anti-symmetrize the 
spin states of these two spin-l/2's. So that anti-symmetrization is done along each 
edge. These anti-symmetrizations ensure that there is no projection on the highest 
bond spin states for every bond. Then we also symmetrize the product of spin-l/2's 
at each vertex (each large circle in Figure [IJ . These symmetrizations preserve the 
correct spin value at each vertex. 

Let us write down the VBS ground state algebraically following this approach. 
We label the particular dot from vertex I connected with some dot from vertex k 
by Ik (correspondingly, that dot from vertex k is labeled by ki). In this way we have 
specified a unique prescription of indices with dots. Then the anti-symmetrization 
results in the singlet state 

|$)fc, = i=(|t)/Ji)fe,-U)/Jt)fe,), (22) 
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Fig. 1. Example of a part of the graph including vertex k with 2^ = 3 and vertex I with 2; = 4. 
Black dots represent spin-1/2 states, which are enclosed by large circles representing vertices and 
symmetrization of the product of spin-l/2's at each vertex. Solid lines represent edges (bonds) 
which anti-symmetrize the pair of connected spin-l/2's. 



where | f) (or i)) denotes spin up (or down) states referring to a basis. The direct 
product of ah these |<i>) singlet states corresponds to all edges (bonds) in our Figure 

m 



(ki) 



(23) 



We still have to complete the symmetrization (circles) at each vertex. We denote 
the symmetrization operator of zi black dots in vertex I by P(0- The action of P(Z) 
on any product of 2; spin-l/2's is 



zi ! terms 



(24) 



where fci, /c2, • . • , fcz, are the zi spin-l/2's (block dots) belonging to vertex I (the 
index also refers to the zi vertices connected to vertex I by an edge) and cr is a 
permutation of the indices: 



1 2 ■■■ zi 
a{l) a{2) ••• a{zi) 



(25) 



All permutations are summed up in (j24p . Then the symmetrization at each vertex 
is carried out by taking the product P(Z) of all vertices. Finally, the unique VBS 
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ground state can be written as 

ivBS)=npwni*)fc'- (26) 

/ (kl) 

Here the first product runs over all vertices and the second over all edges (bonds). 
So that we have constructed a ground state of the AKLT Hamiltonian (PO)) such 
that 

i7|VBS)=0, ^s.+s, (fc, |VBS) =0, V {kl). (27) 

Note that the VBS state in ((26|) is not normalized in general. If the coordination 
number z; is a constant over all vertices in the graph except for boundaries (such 
as in the case of a lattice), then we would have the same spin value at each bulk 
vertex. In that case the basic model is also referred to as the homogeneous model. 

2.4. The Generalized AKLT Model 
2.4.1. The Hamiltonian 

In the generalized AKLT model, relation (|21l) is generalized. We associate a positive 
integer Mki {Mki = Mik) to each edge (kl) of the graph (or each bond of a lattice). 
We shall call Mki multiplicity numbers. Similar to the basic model, the AKLT Hamil- 
tonian describes interactions between nearest neighbors (vertices connected by an 
edge): 

i/-^i?(fc,Z). (28) 

(kl) 

However, the Hamiltonian density is no longer proportional to a single projector as 
in (|19p in general. Instead, it is a linear combination of projectors 

Sk+Si 

H{k,l)^ Cj{k,l)njikJ), (29) 

l<Mki< 2S< , 5< = min{5fc. Si}. 

Projector 7rj(fc, /) is given by p7)) . and Cj{k, Z)'s are arbitrary positive coefficients. 
So that H{k, I) projects the bond spin on the subspace with bond spin value J 
greater than Sk -\- Si — Mki. Physically formation of any bond spin higher than 
Sk + Si — Mki would increase the energy. 

The condition of uniqueness of the ground state was introduced by Kirillov and 
Korepin in^Sl 

2Si=J2Mki, VL (30) 

fe 

Here Si is the spin value at vertex I and we sum over all edges incident to vertex 
I (connected to vertex I), i.e. k are nearest neighbors. The Hamiltonian (|28|) has a 
unique ground state if (|30|) is valid, see S I2.5I . The relation 25"; = zi for the basic 
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model is a special case when all Mki — 1. Also, when M^i = 1, the Hamiltonian 
density (l29t reduces to that of the basic model ([19]). The condition ((30|) can be 
put into an invariant form. Let us define a column vector S, the l^^ component 
of which is associated with vertex I of the graph and equal to Si . The number of 
components is equal to the number of vertices N in the whole graph. Next, we define 
another column vector M with its dimension equal to the number of edges M in 
the graph. The fc*'' and l^^ components of this vector are associated with edge (kl) 
and both equal to Mki ■ The most import ant g eometrical characteristic of the graph 
is the vertex-edge incidence matrix I (see . This is a rectangular matrix with N 
rows and M columns. Each row is associated with the vertex and each column is 
associated with the edge. If the vertex belongs to the edge the corresponding matrix 
element is equal to one, otherwise zero. Then the condition p0[) of uniqueness can 
be re-written as 

2S = i-M. (31) 

For more details we refer to 

2.4.2. The VBS State - Schwinger Boson Representation 

Under condition pOI) or (PT|) . the unique ground state of Hamiltonian (1^51) is referred 
to as the generalized VBS stat e It is c onstructed by introducing the Schwinger bo- 
son representation 1 5 1 25.52,48,76 1 77 1 7 81^ This method uses a pair of bosonic creation 
and annihilation operators (similar to the treatment of the harmonic oscillator prob- 
lem) to realize the SU (2) Lie algebra. 

Define a pair of independent canonical bosonic operators ai and bi for each vertex 
(or site) I: 

[ak,al] = [bk,bl]^5ki (32) 
with all other commutators vanishing: 

[at , ai] = [bk , bi] ^ [uk , k] ^ [ak , bj] ^ 0, \/ k,l. (33) 
Spin operators are represented as 

S+^aX S^^bjai, S[- ^^{ajai-blbi). (34) 

It is straightforward to verify that the SU (2) Lie algebra is realized. To reproduce 
the dimension of the spin- 5; Hilbert space at vertex I, a constraint on the total 
boson occupation number is required: 

Si = ^{ajai+bjb,) = Si. (35) 

i.e. any physical spin state \^)i at vertex I must satisfy Si\4')i = S';|'0)/. In this 
framework the spin state \Si,mi)i such that Sf\Si,mi)i = Si{Si 1)15*;, rn;); and 
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S[\Si,mi)i = mi\Si,mi)i is represented by 

(„USi+mi (uf\Si~mi 

|5,,m,), = yj \;> == |vac),, (36) 

where the vacuum |vac); is annihilated by any of the annihilation operators: 

ai |vac); = bi |vac), = 0. (37) 

As a result, the VBS ground state in the Schwinger representation is constructed 
as 

\YBS) [alb] -blalY''^' \v^c). (38) 

(ki) 

It worth mentioning that this representation shows that for a full graph (each vertex 
is connected to every o ther ver tex by definition) the VBS state coincides with the 
Laughlin wave function 1^1211311 the product runs over all bonds (edges) and 

the vacuum |vac) is the direct product of vacuums of each vertex 

|vac) =(g)|vac)i, (39) 
I 

which is destroyed by any annihilation operators a/ or 6;, V L (Note that [aj. , b]] = 
0, V k,l.) 

To prove that (j38p is the ground state we need only to verify for any vertex I 
and bond (edge) (kl): 

(1) The total power of a] and bj is 2Si, so that we have spin-5; at vertex Z; 

(2) The z-component of the bond spin satisfies 

- ^(E ^'^''^ + E ^'^*'') ^ ^ SI + St < i(E + E Mk'i) 

I'^l k'^k I'^l k'^k 

(40) 

by a binomial expansion. Consequently, the maximum value of the bond spin 
'hi is Mi'k + Efc'/fc A^fc'i)/2 ^ Sk + Si- Mu (from SU{2) invariance, 

see §0 and El). 

Therefore, the state |VBS) defined in ([55]) has spin-S"; at vertex I and no projection 
onto the Jki > Sk + Si — Mki subspace for any bond (edge). As a consequence, 

iJ|VBS)=0, 7rj(fc,/)|VBS) = 0, 
Sk + Si-Mki + l<J <Sk + Si, y {kl). (41) 

The introduction of Schwinger bosons can be used to construct a spin coherent state 
basis (as expected due to the similarity with the harmonic osc illator) in which spins 
at each vertex behave as classical unit vectors, see ii l5.1.2l and I 5 | 25 | 52 | 48 | 76 177j r^j^^ 
coherent state basis converts algebraic computations into classical integrals which 
becomes extremely useful in later sections. 
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2.5. The Uniqueness Condition 



As stated in previous sections, the condition for the existence of a unique VBS 
ground state is 25*; = z; for the basic model and 25"; — J^k ^'^ki '^'-'^ generahzed 
model (the former being a special case of the latter). This uniquen ess c ondition for 
the AKLT model defined on a finite graph or lattice was proved in 

Now let us turn to the proof of the uniqueness condition 25; = Y^^^^ki , i-e. 
the equation 

i?|5')=0 (42) 

with H the AKLT Hamiltonian (pS)) has exactly one solution under the condition 
([30|) or (|3T|) . Note that this expression (j42|) is equivalent to 



TTj{k, /)!*) =0, V (fc, I), Sk + Si-Mki + l<J < Sk + Si (43) 

because of the positive semi-definiteness of every projector ttj and the positive coef- 
ficients Cj. In order to prove the uniqueness condition, we first prove the following 
lemma. 

Lemma 1. All solutions of the equation 

TTjik, OIV-) =0, Sk + Si-Mki + l<J<Sk + Si (44) 
for fixed k and I can be represented in the following form 

|^) = /(at,6t)(45j-a|5t)^^-|vac). (45) 
Here /(a^,5^) is some polynomial in aj,, and aj , bj . 

Proof: For convenience, let us apply the Weyl representation of the SU{2) Lie alge- 
bra. Consider the space of polynomials in pairs of variables xi and yi with coefficients 
in C. We represent operator aj as multiplication on xi and bj as multiplication on 
yi. At site I we have 



d d d d 



2S^ ^xi—-yi — , 2Si^xi—+yi — . (46) 
dxi dyi dxi dyi 

A basis of the (25'; -I- l)-dimensional irreducible representation of spin-S*; can be 
chosen in such a way: 

Vs, = {xf'+™'yf \m = -S,...,S}. (47) 

These monomials with total power 25"/ are clearly eigenvectors of Sf and Si. Now 
let us consider the tensor product of two irreducible representations V5; ® Vs^. . 
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Define the bond spin Jki = Sk + Si, then 



oyk oyi 



d_ _d_ 

dxk ^' dxi 
d d d d 



Jki = yk— + yi—, 



'^Jki = ^kT. — \'Xi- — vkT. — yi-^' 

dxk oxi dyk oyi 

9 d d d , 

2Jki=Xk-^ ha;;- hyfc^ ^yi^- (48) 

dxk dxi dyk oyi 

The tensor product of irreducible representation can be reduced to a direct sum of 
irreducible representations 

S I s 

ys,®Vs,= Vj- (49) 

J=\Sk-Si\ 

Now we construct the highest vector (polynomial) vj of irreducible representation 
Vj with fixed J: 

■^ki'"'J = 0' ^fc;^./ = ■^'"J^ Jkivj = {Sk + Si)vj. (50) 
It must have a total power 2(5'^ + Si), so that the form can be taken as 

— 2^ ^mkm,Xi yi X^ y^ . (OL) 

mk+mi=J 



This form already satisfies the second and third equations in (jSOp . After rearranging 
terms (relabeling indices), the first equation of (|50p becomes 



Jkl'"-J = X! [{^k~mk)C„ikJ-mk+{Sl- J + mk + l)Cmk + i^.J-rak~l] 

yk Vi ■ y^"^) 



nifc=0 

Sfc+rrifc + l Sfc-rrifc-l Si+./-mfc Sj-J+mj, 

'■^k yk -^i yi 

Because of the linear independence of the monomials appearing in ()52p , the coeffi- 
cients must vanish, which yields the following recurrence relation 

Cmk + l,J-mk-l — ~Ti , 1 (-'mk,J-mk- (53) 

bl - J + TUk + I 

The solution to (|53p in terms of Co,j is 

(^\Sk-mk f Sk + Sl- J 

^ ' \ Sk-ruk , 

Cm^J-^m^ = — — — Co, J. (54) 
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Therefore by substituting (j54p into (j5ip and recognizing a binomial expansion, the 
form of vj is found to be 

vj ^ — — — —^k [Xkyi-xiyk) ' 



(-1)' 



Sk + Si-J 
Sk 



(55) 

<xxf''~^'x^'''~'''{xkyi-xiykf', M^Sk + Si-J. (56) 

The over aU constant factor is irrelevant. All other vectors of representation Vj can 
be obtained from the highest vector vj by applications of operator Jjlj . Notice that 
J^i commutes with the factor {xkyi — xiyk) 

[J^i-, Xkyi- xiyk]=0. (57) 

So that all vectors of representation Vj are divisible by 

{xkyi-xiykf, M = Sk + Si-J. (58) 

In other words, any vector (polynomial in Xk, yu and xi, yi) in the vector space 
spanned by Vj has a common factor (|58p . As a consequence, if there is no projection 
on the states with bond spin values 

Sk + Si^ Mki + 1<J <Sk + Si (59) 

after summation of spins Sk and Si {i.e. no projection on Vj with S'^+S'; — Mfej+l < 
J < Sk + Si), then a factor 

{xm - xiykf'^' (60) 

can be isolated, i.e. Any vector in X]j!L|^fc-sf|°' ^ J would have a common factor 
Moreover, this fact is independent of whether we are using the Weyl representation 
or the Schwinger representation of the Lie algebra. Therefore, any solution to (|44l) 
must take the form of (I45p with the factor {a\,b\ — a^^j,)*^''' isolated. Thus we have 
proved Lemma [T] □ 

Now let us use Lemma [1] to prove the uniqueness condition (pO| or ((3T|) . Note 
that (|44t is valid for each bond (fc/), consequently any ground state \^) of the 
Hamiltonian satisfying (l42t and p3|) can be presented in the form 

I*) = YlF{a\b'^){alb\ - at5l)*^-|vac), (61) 

(kl) 

where F{a\b'^) is some polynomial in a''''s and 6'l''s. Now we have to make sure 
that in ((6T|) each vertex (site) should have the correct spin value. By applying 
25*; = {a\ai + b\bi) to the state l^"), we realize that the explicit factor in (|6T1) 
contribute to 25*; (denoting the eigenvalue of 25"/) exactly the value X]{/c;) -^^ki which 
is the sum of powers of a\ and b\ . A comparison with expression ([50]) or (pij) shows 
that if we require this condition 2Si = X](fc/) ^ki, then each site would already have 
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the correct spin value with the presence in (1611) of the exphcit factor only. Therefore 
the degree of the polynomial in variables aj and bj is zero. This is true for every site 
I. Therefore the polynomial F(a^,6^) is a constant which can be removed. So that 
we have proved that the uniqueness condition ((30)) or (1311) guarantees the existence 
(through explicit construction in §[2A2]) and uniqueness of an energy ground state 
- the VBS state. 



3. The Subsystem and Measures of Entanglement 



The VBS states constructed in previous sections ii l2.3l and S I2.4I as ground states 
of AKLT models are highly entangled states. The quantification of entanglement is 
our main subject of study. 

3.1. The Block Density Matrix and the Block Hamiltonian 
3.1.1. The Block Density Matrix and Entropies 

The VBS state (see (|26l) and (|38p ) has non-trivial entanglement properties. The 
density matrix of the VBS state is a projector (a pure state density matrix): 

_ I VBS) (VBS I 



(VBS|VBS) 



(62) 



In order to analyze the entanglement, let us cut the original graph (lattice) into 
two subgraphs (sublattices) B and E. That is, we cut through some edges (bonds) 
such that the resulting graph (or lattice) B U E becomes disconnected (no edge 
between B and E). We may call one of them B, the block, and the other one E 
the environment. The distinction is somewhat arbitrary and the two subsystems are 
equivalent in measuring entanglement. 

Let us focus on the block (subsystem B). It is described by the density matrix 
of the block (obtained by tracing out all degrees of freedom of the environment 
E from the density matrix p ([S^ l: 

PB=tr£;[p]. (63) 

In (j63p and below we use subscript B for block and E for environment. After tracing 
out all degrees of freedom outside the block the density matrix pg is, in general, 
a mixed state density matrix (unless the pure state density matrix p of the whole 
system projects on a product state, which is obviously not our case of the VBS 
state) . Formula (1551) is the definition of the block (subsystem) density matrix and 
it satisfies all three requirements of a density matrix: 

(1) The trace tr^ [ pg ] = 1 and hermiticity = pg follow immediately from those 
of p; 

(2) The positive semi-definiteness is seen by picking up an arbitrary state |'0)-B of 
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the block and realizing that 

b{MPb\^)b = trs [pb\^)b{M] 

= tl B[itVEpMB{lJA] 

= tr[p|V')B(V'|®/ij] >0, (64) 

because of the positive semi-definitencss of p [Ie is the identity of the environ- 
ment). 

The density matrix Pb is a central quantity in description of the subsystem (block). 
It contairL S all correlation functions in the VBS ground state as matrix entries 
| 5 | 46 | 48 | 7 ni (The relation between elements of the density matrix and correlation 
functions is given in ii l5.9l .) It is essential in measuring the entanglement which is 
our main subject. 

The entanglement can be measured or quantified by the von Neumann entropy 



'v N 



-tre [ps Inps ] = - ^ Aln A (65) 



or the Renyi entropy 

^R(a) = -^ln{trB[pS]} = -^ln|VA" , a > 0. (66) 

Here A's are (non-zero) eigenvalues of the density matrix p^. The corresponding 
eigenvector is denoted by |A). i.e. 

Ps|A)-A|A), A^O. (67) 

The Renyi entropy depends on an arbitrary parameter a. If we know the Renyi 
entropy at any a, then we know all eigenvalues of the density matrix. Note that 
the Renyi entropy can be regarded a generalization of the von Neumann entropy 
and reduces to the latter in the limit a — > 1. The von Neumann entropy is a proper 
extension of the Gibbs entropy (in statistical mechanics) and the Shannon entropy 
(in information theory) to the quantum case. (The Shannon entropy measures the 
uncertainty associated with a classical probability distribution. Whereas in quan- 
tum case a density operator replaces a classical distribution.) It was shown by using 
the Schmidt decomposition (Section 2.5 of ^^SJ) lYiat non-zero eigenvalues of the 
density matrix of subsystem B (block) is equal to those of the density matrix of 
subsystem E (environment). So that the two subsystems are equivalent in measur- 
ing entanglement in terms of entanglement entropies, i.e. n[S] = '5'v n[E\ and 
^^[i?] = S'r[£']. This fact has been u sed in obtaining entanglement entropies of 
1-dimensional VBS states as in instead of diagonalizing the density matrix 

directly. We will study the entropies in detail in following sections. 
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3.1.2. The Block Hamiltonian 

The AKLT block density matrix possesses certain characteristic properties which 
distinguish the VBS states from others. We shah show in §[231 that the spectrum 
of the density matrix contains a lot of zero eigenvalues. In order to understand 
this and give the subsystem (block) a more complete description, we first introduce 
the Hamiltonian of the subsystem (called the block Hamiltonian). 

The block Hamiltonian is the sum of Hamiltonian densities -ff (fc, /) with 
both k £ B and I € B, i.e. nearest neighbor interactions (bond terms) within the 
block B: 

Hb= ^(^'0: keB, leB. (68) 

{kl)eB 

Here H{k, I) is given in (|19p for the basic model and (I^Hl) for the generalized model, 
for k and I connected by an edge (bond). In ()68|) no cut edges are present (boundary 
edges between subgraphs B and E removed). In other words, the block Hamiltonian 
is the internal interactions of the subsystem B. This Hamiltonian has degenerate 
ground states because uniqueness conditions ([2T|) and (pO]) are no longer valid. 

Let us discuss the degeneracy of ground states of the block Hamiltonian (ISS]) . 
Let us denote by NgB the number of vertices on the boundary dB of the block B. 
The boundary consists of vertices (sites) with several incident edges (bonds) being 
cut, see Figured] The degeneracy of ground states of Hb (which is abbreviated by 




Fig. 2. Example of the cutting for the basic model. The curved double line represents the bound- 
ary between the two subgraphs. We have the block B on the left and the environment E on the 
right. Solid lines represent edges (bonds) while dashed lines represent cut edges (cut bonds). Each 
dashed line connects two dots. All vertices in the figure belong to the boundary of B or i? because 
of the presence of one or more cut incident edges (dashed lines). 
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deg) is given by the Katsura's formula 



deg = Y[ 



^ Mki + 1 

iGdB L \k<£dE ) 



{kl) G {cut edges}. (69) 



Here dB denotes vertices (sites) on the boundary of the block B and dE are vertices 
(sites) on the boundary of the environment E. In (1691) we have Nqb terms in the 
product. Formula (15^ is valid for both the basic and the generalized model. For 
the basic model all M^i — 1, including those corresponding to cut edges. Take, for 
example, a particularly simple case that each vertex on the boundary of the block 
dB was connected to exactly one vertex on the boundary of the environment dE. 
Then the degeneracy deg — 2^°^. A general proof of formula is given in the 
next section ? I3.2I . The subspace spanned by the degenerate ground states is called 
the ground space, with the dimension given by deg in (|69p. We emphasize at this 
point that the block B should contain more than one vertices, otherwise we have 
a trivial case that the block Hamiltonian vanishes Hb — and the whole Hilbert 
space become the ground space. We discuss the density matrix for a single vertex 
block at the end of §[33]. The spectrum of the density matrix pg is closely related 
to the block Hamiltonian. The density matrix is a projector onto the ground space 
multiplied by another matrix. We shall prove this statement for an arbitrary graph 
or lattice in ^\3.3\. 



3.2. The Degeneracy of Ground States of the Block Hamiltonian 

We prove Katsura's formula (IMl) for the degeneracy of ground states of the block 
Hamiltonian. The proof applies to both the basic and the generalized models. The 
block Hamiltonian is defined in Let us first look at the uniqueness condition 
([50]) . (The condition (pij) for the basic model is a special case of this general one.) 
For an arbitrary vertex (site) / in the block B, the condition can be written as 

2Si = Mki = ^'^ki + ^ Mfci, / e B. (70) 

k keB keOE 

Note that the sum over vertices k e dE is outside the block B. These terms are 
only present for boundary vertices I e dB. Expression (j70p is valid for any vertex 
in the block (for a bulk vertex the last summation vanishes). Next we define the 
block VBS state 

|VBSArJ= n {^bl-biaiy'"' \vac), keB, leB. (71) 

Here edge (bond) {kl) lies completely inside the block B. Now an arbitrary ground 
state of the block Hamiltonian Hb takes the following form (see Lemma [1] in §[13] 
which explains the appearance of the factor {a\.b\ — b\,a[Y^''^ in (|7l]) above): 

~Nqs terms 

iG)= n f^^iA) 

I e dB 



iVBSjvJ, (72) 
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where f{aj, bj) is a polynomial (it may depend on the vertex I) in aj and bj and the 
product runs over all boundary vertices (with the total number denoted by Nqb)- 
The degree of this polynomial is equal to J^kedE -^^ki- (Each term in the polynomial 
must have the same total power J^kedE ^^ki of aj and bj.) It is straightforward to 
verify that |G) in ((72)) is a ground state: 

(1) The power of aj and bj in IVBS^Vb) is J2keB (see ((7T|) ). so that the total 
power of aj and bj in ([71]) is J2keB ^ki + J2kGdE ^ki = '^Si according to ([TO]). 
Therefore, we have the correct power 25*; of the bosonic operators aj and bj for 
each vertex / in the block B (constraint (PSI) is satisfied); 

(2) There is no projection on any bond (edge) spin value greater than or equal to 
Sk + Si — Mki + 1 because of the construction of the block VBS state (|7T|). (One 
could use the same reasoning as in §[211]). 

Therefore the degeneracy deg of the ground states of Hb is equal to the number 
of linearly independent states of the form (17^ . Since a|'s and bj's are bosonic and 
commute, the number of linearly independent polynomials f{aj, a\) for an arbitrary 
I is equal to its degree plus one, i.e. 



kedE ) 



_ 1, yiedB. (73) 

\ke ' 

So that the total number of linearly independent polynomials of the form 

Nob terms 

n /K''^?) (74) 

; e OB 

is the product of these numbers ([75]) for each I e dB. Finally, the ground state 
degeneracy of the block Hamiltonian Hb is (Katsura's formula) 



deg= Y[ 
ledB 



kedE / 



1 



(75) 



^feec 

In the case of the basic model all M^i = 1, formula ()75p has a graphical illustration, 
see Figure [2] We count the number # of all cut edges (dashed lines) incident to one 
boundary vertex of the block, then add one to the number The degeneracy is 
the product of these (# + l)'s for each boundary vertex. 



3.3. General Properties of the Density Matrix 

The reduced density matrix from a VBS state has important and special spec- 
trum structures which are universal among AKLT models. Let us denote by Nb the 
number of vertices in the block B. Then the dimension dim of the Hilbert space of 
the block B is equal to n/(2'S'/ + 1) with / G B, which is also the dimension of the 
density matrix p^. The value is 

dim = [z, + 1] , (76) 
leB 
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for the basic model and 

dim = 

leB 

for the generahzed model. In both expressions ([75)) and (|77l) we have Nb factors 
in the product. Take, for example, a particularly simple basic model in which each 
vertex is connected with the same number z of vertices, including those correspond- 
ing to boundary vertices (such as in the case of a lattice). Then the dimension 
dim = (z + 1)^^. The density matrix pg would have dim number of eigenvalues. 
However, most of the eigenvalues are vanishing and is a projector onto a much 
smaller subspace multiplied by another matrix. To prove this statement, wc define 
a support to be the subspace of the Hilbert space of the block B with non-zero 
eigenvalues, i.e. it is spanned by eigenstates of pg with non-zero eigenvalues. The 
dimension of the support is denoted by D. Then We have the following theorem on 
the structure of the density matrix pg (Assuming that the block have more than 
one vertices, i.e. Ng > 2, so that Hg is not equal to zero identically): 

Theorem 2. The support of pg Ii63\) is a subspace of the ground space of the block 
Hamiltonian Hg i68\) . 

Proof: To prove the theorem, we recall that H = J2{ki)eB ^(^^ each H{k, I) 
is a sum of projectors (P^ . We have H{k, I) > 0. Then the construction of the VBS 
ground state P5|) and guarantees that there is no projection onto the subspace 
with higher bond spins {J > Sk + Si — Mki + 1) for any bond (edge). (See ? I2.5I for 
the proof.) Therefore, 

H{k,l)\VBS) =0, y (kl). (78) 

In particular, this is true for bonds (edges) inside the block B, i.e. both fc G _B and 
I £ B. Now, from the definition of pg in (|63l) . we have 

H{kJ)pg= H{k,l)tTE[p] 

_ g(fc,Otrg[|VBS)(VBS|] 
~ (VBSjVBS) 
_ trg [g(fc,0|VBS)(VBS|] 
(VBS|VBS) 

= 0, keB, le B. (79) 

In the last step of ([79]) we have used ([78l) and the fact that bond (edge) (kl) lies 
completely inside the block B so that H(k, I) commutes with the tracing operation 
in the environment E. Equation (|79|) is true for any bond (edge) in _B, so that 



Ml, 



(77) 



HBPg= ^ H{k,l)pg=0, keB, l£B. 

{kl)eB 



(80) 
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If we diagonalize the density matrix 

where |A) is the eigenstate corresponding to eigenvalue A. Then ([M]) can be re- 
written as 

HbJ2\\X){X\=Y.XHb\X}{X\=0, (82) 

Note that {|A)} is a hnearly independent set. Therefore the solution of (|82l) means 
that 

Hb\X)^0, A^O. (83) 

Expression (j83l) states that any eigenstate of (with non-zero eigenvalue) is a 
ground state of Hb- As a result, we have proved Theorem [2] that the support of Pb 
is a subspace of the ground space of Hb, so that D < deg. The density matrix takes 
the form of a projector multiplied by another matrix (a constant matrix depending 
on non-vanishing eigenvalues) and the projector projects on the ground space. □ 

Also, it is clear from expressions (|69|) and (|76|) . (f77| that deg < dim (dB C B 
so that Nqb < Nb)- Usually, deg is much smaller than dim because the former 
involves only contributions from boundary vertices (sites) of the block while the 
latter also involves contributions from all bulk vertices (sites). Then as a corollary 
of Theorem [21 we have 

D <deg < dim. (84) 

If the block B consists of only one vertex with a spin-S*, then we conjecture that 
it is in the maximally entangled state. The support has dimension D = 2S + 1. 

4. The One dimensional Spin— 1 Model 

One of the most simple models is defined on a 1-dimensional lattice with spin-l's 
in the bulk and spin-l/2's at both ends. We shall denote by Sj the vector spin 
operator at site j {j = 0,1, . . . , N + 1). The Hamiltonian is 

1^-1/ -| 2\ 

^=2^[^^' + 3 ■ + 3 j + (0, 1) + ^3 (iV,iV + 1). 

j — 1 

(85) 

Each bulk term is a projector 7r2 onto the states with bond spin-2. The boundary 
terms describe interactions of a spin-1/2 on the boundary and a spin-1 in the 
bulk. Each term is a projector onto the states with bond spin-3/2: 

Tr.{0,l)^l{l + So-Si), Tr3{N,N + l) = l{l + SN-SN+i)- (86) 
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The choice of boundary terms guarantees the uniqueness of the ground state. As 
mentioned before, if we have spin-1 at every site in (|85p instead, the ground state 
would become 4-fold degenerate. 

In this section we study the entanglement of the unique VBS ground state of 
this 1-dimensional spin-1 model. As to be shown below, the density matrix of a 
block of L contiguous spins is diagonalizable. It has four non-zero eigenvalues: 



A„ 



i(l-f 3(-i)^), a = 0; 



a = 1,2,3. 



(87) 



These eigenvalues depend on the length L of the block subsystem and are indepen- 
dent of the size of the whole spin chain. The von Neumann entropy and the Renyi 
entropy derived from these eigenvalues are 



ln4-i(l + 3(-l)^)ln(l + 3(-l)^) 
-J(l-(-^)^)ln(l-(-J)^) 



1 — a 



■In 



J(l + 3(-i)^) 



r(i-(-o)'') 



Note that the parameter a in the Renyi entropy should not be confused with the 
label a for the eigenvalues in ([87|) . 



4.1. The VBS Ground State 

Given the Hamiltonian (|85l) . we are going to use the pictorial method (see ^^ 12.3.21 ) 
to construct the unique VBS ground state. In order to represent the state, we first 
introduce the following notation for convenience '^Sl 

|a) = (-l)i+'^°'«/(8)ac|0), a = 0,1, 2, 3 (89) 

where ao = I (2-dimensional identity), aa=i,2,3 are standard Pauli matrices and 
|0) = —(I ti) ~ I it))/'s/2 is the singlet state (antisymmetric projection) of two 
spin-l/2's. (It corresponds to the antisymmetrized state |<I>) in ii l2.3.2l .') These four 
states (|89p (called maximally entangled states) form an orthonormal basis of the 
Hilbert space of two spin- 1/2 operators. 

The spin-1 state at each site is represented by a symmetric projection of two spin- 
1/2 states given by ([Ml) for a = 1, 2, 3. Let us take the j**^ site for example, see Figure 
m The two spin-l/2's are labeled by {j,j) (from left to right, respectively). Then the 
spin-1 states are prepared by projecting these two spin-l/2's (4-dimensional space) 
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11 j j 7 + 1 J+l N N N + l 



Fig. 3. Graphical representation of the VSB ground state for the 1-dimensional spin-1 model: 
Each spin-1 is realized by a pair of spin-l/2's which are represented by small black dots in the figure. 
The pair of spin- 1/2 states at site j are labeled j, j. The solid lines connecting two neighboring 
dots (j and j + 1) represent anti-symmetrization of two spin-l/2's; The large circles enclosing two 
dots (j and j) represent symmetrization at each site. The boundary spin-l/2's are labeled and 
A'' -f 1 in consistency with the prescription. 

onto a symmetric 3-dimensional subspace spanned by 



l% = ^(lt>,lt>j-U>,U)j), 

|2),j = ^(|t>,|t>j + U>,U>5), 



V2' 



|3),j = -=(|t>,U>j + U>,|t>5). (90) 



Thus the two ending spin-l/2's are labeled as site and iV -|- l, co nsistently (Figure 
[3]) . The unique VBS ground state in this representation is I I I ° l 

N 

|VBS) - ( (g)P,j I |0)oi|0)i2---|0)^^+i. (91) 

Here Pjj projects two spin-1/2 states onto a symmetric subspace, which describes 
spin-1. Using basis (15^1) . we have 



^|a),j(a|. (92) 



This projector P^j serves the same purpose as the symmetrization operator P{1) 
in ? I2.3.2I and their results acting on a product state of spin-l/2's only differ by a 
normalization. We use the projector P^j for convenience here. 

A crucial step (see is that the ground state can be expressed in a 
different form using the following identity 



-1 ^ 

\0)ab\0)bc = — E ^ (^a)c] \0)ac (93) 



for arbitrary labels (indices) A, B and C. This identity (1^ can be verified by direct 
calculation and comparison. Repeatedly using relation (|93|) . the product of |0)'s in 
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(Plj) can be re-written as 

|0)oi|0)i2---|0)^^r+i (94) 

N 3 

\ai) ■ ■ ■ law) [Jo «) (cra„ • • -CTaJAT+i] |0)o7V+l- 

ai ^ctN—O 

Then by projecting onto the symmetric subspace spanned by \a — 1,2,3), the 
ground VBS state ^ takes the formal! 

1 ^ 

= ^ Iai)---|"A') [4® (f^a^-'-^aJAT + i] |0)oN+l- (95) 

Ql ■ ■ ,QiV — 1 

Note that this ground state (|95|) is normahzed and we have re-written the overah 
phase for it has no physical content. 

4.2. The Block Density Matrix 

Given the ground state in the form ()95p . we obtain the density matrix of a block 
of L contiguous bulk spins starting at site k by tracing out spin degrees of freedom 
outside the block using basis (|89|) : 

PL = t^-o,i,...,k-i,k+L,...,N,N+i [ I VBS) (VBS| ] . (96) 

(Note that we use subscript L to emphasize the dependence of the density matrix 
on the size of the block instead of using the general B as representing 'block'.) In 
taking the partial trace, we will encounter the following expression in calculation 

3 

^ /®(tTa„---^Tai)|0)(0|/(8)(cr„, •••a„J, n>l (97) 

given Iq — |0)(0|. To solve this (|97l) . we introduce iterative coefficients An and i?„ 
and write 

3 

/„=A„|0)(0KB„^|/3)(^|, n>l. (98) 

13=1 

Then from (I97[) we could write down the expression of /„+i in terms of An and Bn- 
Comparison of coefficients yields the following iteration relation 

An+l = Wn, Bn+1 = A„ + 2B„, n > 1 (99) 

with Aq = 1 and Bq = 0. The solution to dM]) is 

A„ = J(3" + 3(-l)"), B„ = i(3"-(-l)"). (100) 
As a result, we have found that 

In^l (3" + 3(-l)") |0)(0| + \{T~ (-1)") ^ |/3)(/3|. (101) 

/3=1 
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Using (|10ip . it is straightforward to take the partial trace in (|M|) . The result is 
independent of the starting site k and the total length TV (see^l^. (So that the den- 
sity matrix is translational invariant though the whole spin chain Hamiltonian does 
not have complete translational invariance because of the boundary conditions.) We 
choose A; = 1 {i.e. re-label the indices of sites for notational convenience) so that 
the density matrix reads 

1 ^ 



3^ 

Q.a — 1 



(102) 



4.3. Ground States of the Block Hamiltonian 



The block in 1-dimension is L contiguous bulk spins. The block Hamiltonian Hb 
by definition (|68p reads 

= ^ E (^^ ■ ^^+^ + \ (^j- ■ + 1) ■ ^^^^^ 

Any linear combination of states of the following form 

|G;xi,Xl)^ (^^'') '^^^i'°^i2'"^23---|0)r3TLlXL)L (104) 

is a ground state of the block Hamiltonian (|103p . In (|104p we have made notation 
|x) = I t or i) represents the two spin-1/2 states and Pjj is defined in ([92]). Let us 
make a particular linear combination of these |G; xi, Xl) states using (|89p and write 
the four (a = 0, 1, 2, 3) linearly independent ground states of the block Hamiltonian 
(fT03l) as follows 




|VBS;a)^ | (X)P,j | |a)£i|0)i2|0)23 • • • |0)^^. (105) 

Note that we have changed the label G to VBS and these 4 states in (|105l) are called 
degenerate VBS states. Now we go through the same steps as from (|9T|) to (j95|) . the 
resultant form of the four ground states (a = 0, 1, 2, 3) is 

3 

|VBS;a)= ^ | ai) ■ ■■ \aL) (ailo-Q (Ki (o-q^_i • • • ctq J |0). (106) 

ai,--- :aL = l 

Again we have re-written the overall phase for simplicity. These four states are 
orthogonal, and the normalization is given by (the calculation is similar to that of 
A„ and B„ in (|M|) ) 

i(3^ + 3(-l)^), a = 0; 
(VBS;a|VBS;a) = <( (107) 
K3^-(-l)^), a = 1,2,3. 
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4.4. Spectrum of the Density Matrix 

According to Theorem [21 the eigenvectors corresponding to non-zero eigenvalues 
of the density matrix (|102l) are degenerate ground states of the block Hamiltonian 
(I1G3|) . These are exactly the degenerate VBS states found in (|106p . Let us apply 
to |VBS;a) and use orthogonality of the \a) states. Then we obtain 

1 ^ 

Pi|VBS;a> = ^ J2 • • • I^l) C„,...„, (108) 

Ql,--- ,Qi = l 

with coefficient 

3 

Cci-ai, = ("Lka® (fTa'^_^ •••cr„;)|0) (109) 

•(0|/ (g) (ct„'^ • • • a„^)/ (cr^^ • • • CT„ J|0). 

Using the same method of induction as in obtaining An and i?„ in (|M| . we have 

3 3 

with 

[- 1(3^-1 +3(„l)L-i), /3 = 0; 
Ap ^ \ (111) 
I i(3^-i- (-1)^-1), /3 = 1,2,3. 

Therefore the coefficient Cai - UL defined in (|109p can be simplified as 

3 

)|o). 

(112) 

Straightforward calculation using multiplication rules of Pauli matrices shows that 
pi2p can be further simplified as 

Safi{aL\I (g {(JaL-i ■ ■ ■ <^ai)\0) (113) 
3 

+ (Ao + 2Ai)(l - <5„,o)(W (0| - i ^ )|o) 

where CaaLP is the totally antisymmetric tensor of three indices with 6123 = 1. By 
realizing that 

3 

5aaL{0\ - i^ eaai./3(/3| = (OlcTa^CTa (ailcTa ® I, (114) 

/3 = 1 

we have reached the final form of the coefficient Cq,j...q,^ such that 

C„i...„^ = [3Ai(5„,o + {Ao + 2Ai){\ - So^^o)] {aL\(Ta ® {(^ar^-i •••c^ai)|0). (115) 
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As a result, we plug (IllSp into (|108p and find that 

IVRQ \ 3^i<5^^o + (^ + 2v4i)(l- J^,o) 
Pi|VBS;a) = (116) 

3 

^ |ai) • • • Itti) (ailcTa «) (tTai,-! • • -CTaJIO). 

Ql,--- ,CtL = l 

By comparing with (|106p . we find that (|116p is exactly the statement that |VBS; a) 
{a = 0, 1, 2, 3) are eigenvectors of the density matrix p^: 

Pi|VBS;a) = Aa|VBS;a), a = 0,1,2,3 (117) 

with eigenvalues 

. 3^1,5^.0 + (Ao + 2Ai)(l-,5„,o) J 3(l + 3("l)'^)> « = «; 
Aq - - < (118j 

I i(l-(-i)^), a = 1,2,3. 

These numbers obtained in (jllSp are exactly the eigenvalues found in fQj. 
spin-1, and are consistent with our later explicit expression for eigenvalues in the 
more general case, see (|20ip in ii l5.6l . 

We can also prove explicitly that any other eigenvectors of orthogonal to the 
set {|VBS; a)} have zero eigenvalue. Note that a complete basis of the Hilbert space 
H of the block of spins can be chosen as 

{|ai)---|aL)}, a = 1,2,3. (119) 

The subspace H4 with non-zero eigenvalues is panned by {|VBS;a)}, as we have 
already shown. The Hilbert space can be reduced into a direct sum 

H = H4®H<[,. (120) 

We will show that the subspace H<[, orthogonal to H4 is a subspace of vanish- 
ing eigenvalues. Mathematically, this means that for an arbitrary basis vector 
• • • we shall have 

Pi(I-P4)|/3i)---|/3L)=0, (121) 

where I is the identity of H and P4 is the projector onto H4: 

3 

1= X! • • • |aL)(ai| • • • (ail, 

Ql,--- ,Ql,=l 

_^ |VBS;a)(VBS;«| 
hi (VBS;a|VBS;a)- 

By taking expressions (|102l) . (|122l) . (|117l) . and realizing that 

t (VBS;a?VBS;a) = g ""^ = ^ ^ ^' ^'''^ 
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we find the left hand side of (|121l) being equal to 

Pi(I-P4)|/3i)---|/3L) (124) 
1 ^ 

= ^ H |ai)---|aL)(O|[/®(a0, •••a^J, /®(aa^---a„J]|0). 

We use multiplication rules of Pauli matrices to write the two terms within the 
commutator in (|124p as 

/ {ap, • • • a/3 J = e"^'-^h ®ap, /? = 0, 1, 2, 3; 

/(g)(cr„^ =e'^(")/®cr„, a = 0,1,2,3. (125) 

Here e'^^^-* and e'^*^"-' are two phase factors. Then the commutator is 

[ / ® {ufi, ■■■on,), I® (a„, • • • a„ J ] = ® [ , a„ ] . (126) 

There are two possibilities: 

(1) a — or at least one of the two is equal to zero, then and CTq. commutes; 

(2) a 7^ /? 7^ 0, then [cr^ , CTq ] — 21X^7=1 ^Pa-^'^ii but we still have (0|/ ® <y^\S) = 
(0|7) = 0. 

Therefore, the factor (0| \l®{ap^ ■ ■ ■ tr/j^) , I®{aai^ ■ ■ ■ CqJ ] |0) in ()124p is identically 
zero. So that we have proved ()12ip . Therefore H$ is a subspace with only zero 
eigenvalues. 

4.5. The Large Block Limit 

It is interesting to study the large block limit that i — )■ oo. We recognized from 
(|118p that all four eigenvalues approach the same limit 

A, = i, L ^ (X3. (127) 

As a result, the von Neumann entropy coincides with the Renyi entropy in the 
numerical value and both equal to 

S-^ N = S'R(a) = ln4, L -> oo. (128) 

The limiting density matrix (limL^oo Pl = Poo) proportional to the projector 
P4 (|122p which projects on the 4-degenerate ground states (the ground space) of 
the block Hamiltonian, i.e. 



(129) 
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5. The One dimensional Spin-S Homogeneous Model 

In l-dimension, if all bulk spins take the same integer value S, the model is called 
the homogeneous model. The system consists of a linear chain of N spin-S"s in the 
bulk, and two spin-5'/2's on the boundaries. Let Sj denotes the vector spin operator 
at site j (j = 0, 1, . . . , TV + 1). The Hamiltonian is 

N-l 2S 

E CjU,J + 1)7Tj{j,j + 1) + H{0,1) + H{N,N + 1), (130) 

j=l J=S+1 

where the projector 7rj(j, j + 1) projects the bond spin Jj,j+i = Sj + Sj+i onto the 
subspace with total spin J {J = S + 1, . . . , 2S). Physically formation of bond spins 
with these values would increase energy. The boundary terms describe interactions 
between a spin-S'/2 and a spin-S": 

3S/2 

F(0,1)= ^^./(o,i)^,/(o,i), 

J=S/2+l 
3S/2 

H{N,N + l)= Cj{N,N + l)7Tj{N,N + 1). (131) 

J=S/2+l 

Coefficients Cj{j,j + 1) can take arbitrary positive values. This model is a special 
case of the generalized model in l-dimension with all multiplicity number Mjj^i — 
S. 

We study the entanglement of the unique VBS ground state of the Hamiltonian 
(|130[) in this section. The density matrix of a block of spins and the entropi es were 
calculated by H. Katsura, T. Hirano, Y. Hatsugai and their collaborators inSSESl 

5.1. The VBS Ground State 

5.1.1. The Construction of the VBS State 

According to the general approach in i; l2.4.2l . the unique VBS ground state of the 
Hamiltonian (|130p is constructed in the Schwinger representation as ^ 

N g 

IVES) ^ n l^ac), (132) 

where a\ are bosonic creation operators and |vac) is destroyed by any of the 
annihilation operators a, b. Recall that these operators satisfy [a^ , ] = [6^ , fej ] = 
6ij with all other commutators vanishing. The spin operators are represented as 
Sj' = Ojbj, S~ = &]aj, Sj = {a^jaj — b^jbj)/2. To reproduce the dimension of the spin- 
s' Hilbert space at each site, an additional constraint on the total boson occupation 
number is required, namely (a]aj + b'jbj)/2 — S. More details and properties of the 
VBS state in the Schwinger representation can be found in S I2.4.2I and l ^ l ^ l ^ . The 
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pure state density matrix of the VBS ground state (|132[) is 

|VBS)(VBS| 



(VBS|VBS) 



(133) 



We will discuss the normalization (VBS | VBS) of the VBS state after introducing 
the coherent state basis. 



5.1.2. The Coherent State Basis 

In order to calculate the normalization of the VBS state (I132p and later the density 
matrix of the block, it is convenient to introduce a spin coherent state basis. As 
shown in^, this basis represents the quantum spins in the model in terms of classical 
unit vectors. We first introduce spinor coordinates 

(u,i;) = ^cos^e'^,sin^e"'^^ , < 6* < tt, < < 27r. (134) 

Then for a point 17 = (sin 9 cos 0, sin 6 sin (f>, cos 9) on the unit sphere, the spin-S* 
coherent state is defined as 



ua 



2S 



\n) = ^ ' Ivac). (135) 

Here we have fixed the overall phase (a J7(l) gauge degree of freedom) since it has no 
physical content. Note that (|135l) is covariant under SU{2) transforrns (see ? l5.8n . 
The set of coherent states is complete (but not orthogonal) such that 



25 + 1 
47r 



s 

dn\n){n\= J2 \S,m)(S,m\ = hs+i, (136) 

m=-S 

where \S,m) denote the eigenstate of and Sz, and hs+i is the identity of the 
(25 + l)-dimensional Hilbert space for spin-5. To prove ()136|) . we expand the ex- 
pression ()135|) (see also ([Ml)) 



m— — 5 V ^ 

Then, by substituting (|137l) into (|136|) and realizing that 



J dClu 



S-\-m^S-m^*S+m'*S-m' _ {S + m)\{S — m)\ 



V U V = (2S +1)1 ^TTdmm', (138) 



the completeness relation (jl36p is then established. This relation ()136p can be used 
in taking trace of an arbitrary operator. 
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5.1.3. Normalization of the VBS State 

The VBS state |VBS) defined in ()132|) is not normalized. Using the coherent state 
formalism (jl35p and the completeness relation (|136p . we express the norm square 
as 



(VBS|VBS) 



(139) 



{s + iy. 



1 2 



An 



(2^+1 )! 
An 



1 N 



N+l \ N r 

\ 3=0 / J=0 L 



n S 



where we have used 



-,3+1, S-l 



(140) 



In order to carry out the integral in (|139p , we consider the expansion of the function 
[^(1 — a;)] in terms of Legendre polynomials 



(141) 



as follows 



^(l-x) 



(142) 



(=0 



The coefficient Ci is derived by using the orthogonality of Pi and repeatedly inte- 
grating by parts 

21 + 1 "'^ 



C, 



2 

2/ + 1 



dxPiix) 



-1 
1 



i(l-) 



1 



^^2^ U?) 



2(1--) 



(2^ + 1)5' 



2s+i+ni{s - ly. 

{-iy{2i + i)s\ 
2s+i+H\{s-iy 

(-1)'(2; + 1)5'!5'! 

{s - iy.{s + 1 + ly 



dx{x^ -lyii-x) 



S-l 



da;(l-2;)^(l + a;)' 



(143) 



Having expansion coefficients (|143p and by replacing x with (ij ■ flj+i, the fac- 

tor i(l — Ctj ■ f^j+i) under the integral of (I139P can be expanded in terms of 
Legendre polynomials and further in terms of spherical harmonics by further using 



An 



(144) 
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The final result is 



n S 



-(l-^-%+l) 



1=0 



s + 



1=0 



=-i 



with coefficients X{l,S) given by 

AG, 5): 



(~1)'5!(^+1)! 

(s - ly.is + 1 + ly: 



(145) 



(146) 



Now we expand ^(1 — Clj ■ (ij+i) 



in terms of spherical harmonics as in (jl45l) . 

then integrate from (1q to Ct^+i- We notice by using the orthogonality of spherical 
harmonics that each integral contributes a factor of 47r/(5+ 1) except the last one. 
For example, 

/■ „ ri . . 



47r 



S- 



-^A(/,5) ^ VI^YC^iCli) / cli]oi^i™(f2o)ro*o(f^o) 



An 



in Y*oin^) 



An 



(147) 



5'+r ^"^^^ 5+1- 

The last integral over ^n+i contributes simply a factor of An. Consequently, the 
norm square (|139p is equal to 



(VBS|VBS) = 



(2S_+1)\ 
5"+ 1 



S'!(S'+ 1)!. 



(148) 



5.2. The Block Density Matrix 

We take a block of i contiguous bulk spins as a subsystem. Now we calculate 
the block density matrix in the VBS state (|132p . By definition, this is achieved 
by taking the pure state density matrix (I133P and tracing out all spin degrees of 
freedom outside the block: 

Pi = tro,i,...,fc-i,fe+L,...,Ar,Ar+i [p] , l<fc, k + L - 1 < N. (149) 

Here the block of length L starts from site k and ends at site k -\- L — 1. pj^ is no 
longer a pure state density matrix because of entanglement of the bloc k with the 
environment (sites outside the block of the spin chain). It was shown inSHESIthat 
entries of the density matrix are multi-point correlation functions in the ground 
state. We give the proof of this statement for our spin-5 case in ^ \5.9\ . 



June 11, 2010 0:33 WSPC/INSTRUCTION FILE AKLT 



36 V. E. Korepin & Y. Xu 



Using the coherent state basis (|135p and completeness relation (|136p , can be 
written asSH 



Pl = 



(150) 



/ 


k-l N+l 

n n d^^. 

_j=0 j=k+L 


fe-2 N 

n n 

j=0 j = fe+L 


i(l-%.f2,+i) 


s 

B^YBSl){YBSl\B 




■(25+ I)!' 


L 

J 


r 


N+l 

J=o 


N 

n 

3=0 




s 



Here the boundary operator B and block VBS state |VBSl) are defined as 

s s 
B = (ufe-i&fe - Vk-iak) {ak+L-iVk+L - bk+L-iUk+L) , 

fe+L-2 g 

iVBSi)^ n - Ivac), 

j=fe 



(151) 
(152) 



respectively. Note that both B and |VBSl) arc SU{2) covariant (see ? l5.8n . The 
expression p50p can be simplified. We can perform the integrals over U,j [j = 
0, 1, . . . , fc — 2, fc + L + 1, . . . , A^, + 1) in the numerator and all integrals in the 
denominator (see § 15.1.31 ). After integrating over these variables, the density matrix 
Pi^ turns out to be independent of both the starting si te k and the total length 
N of the spin chain. This property has been proved in ^^^l for spin 5=1 (using 
a different representation, namely the maximally entangled states, see S I4.2I ) and 
generalized in ^8 for generic spin-S*. Therefore, we can choose fc = 1 (a relabeling 
for convenience) and the density matrix takes the form 



Pl = 



S + l 



with 



(25+1) 



B 



{S + l) 
(47r)2 



dl^od^L+i Bt|VBSL)(VBSL|B 



IVBSl) 



L-l 



n ( a^^'' 



(153) 

(154) 
(155) 



The state |VBSl) is called the block VBS state. The last two integral of (|153l) can 
be performed, but we keep its present form for later use. 



5.3. Ground States of the Block Hamiltonian 
5.3.1. Degenerate Ground States 

In order to describe the eigenvectors and spectrum of the density matrix (|153p . 
we first study the zero-energy ground states of the block Hamiltonian. The block 



June 11, 2010 0:33 WSPC/INSTRUCTION FILE AKLT 



Entanglement in Valence- Bond- Solid States 37 



Hamiltonian is a collection of interacting terms within the block, i. e. 

L-l 2S 

Hb = J2 J2 C,jTTjij,j + l). (156) 

j=l J=S+1 

Now we define a set of 5 + 1 operators covariant under SU{2) (see ^ IS.SI ) 

= (uaj +u6iy'' (a|fe[ - &la[)^ (w4 ^ < J < S". (157) 

These operators act on the direct product of Hilbert spaces of spins at site 1 and 
site L. Then the set of ground states of (|156p can be chosen as 

\G;J,n)=Al\YBSL), J = 0,...,S. (158) 

Any state \G;J,Cl) of this set for fixed J and is a zero-energy ground state of 
(|156p . To prove this we need only to verify: 

(1) The total power of a\ and b\ is 25*, so that we have spin-S at the first site; 

(2) The bond spin value satisfies —S< Jf 2 = 5*^ + 51 < S* by a binomial expansion, 
so that the maximum value of the bond spin J1.2 is S (from SU{2) invariance, 
see §[5l] and IS). 

These properties are true for any other site j and bond (j, j + respectively. There- 
fore, the state |G; J,Cl) defined in (|158p has spin-S at each site and no projection 
onto the Jj,j+i > S subspace for any bond. 



5.3.2. Degenerate VBS States 

The set of states {\G;J,Ct)} depend on a discrete parameter J as well as a continuous 
unit vector Cl. States with the same J value are not mutually orthogonal. It is 
possible also to introduce an orthogonal basis in description of the degenerate zero- 
energy ground states. This new basis could be used in determining the rank and 
the completeness of the set {\G; J, tt)} . For notational convenience, we define 

= V(J + Af)!(J-M)!' = ^iS + my.{S-my: ^'^^^ 

These two variables transform conjugately with respect to one another under SU{2). 
(See S I5.8I for more details of transformation properties.) Variable Xjm has the 
following orthogonality relation 

47r 

dn x*j^tXjM' = (^7^ ^MM' . (160) 

Operator ■0sm ^ spin state creation operator such that 

4„|vac) = |^,m). (161) 
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With the introduction of these variables (|159p , the operator A j defined in (|157p can 
be expanded as (see Chapter 9 



J mi —M 

M— — J mi,mL 

where (S'/2, mi; 5/2, J, M) are the Clebsch-Gordan coefficients. Note that 
mi ^'^'^ defined in the Hilbert spaces of spins at site 1 and site 

L, respectively. We realize that the particular form of the sum over mi and in 
(|162p can be identified as a single spin state creation operator 

mi+TTii— Af 

*JM= E (^/2,™i;5/2,m2|J,M)4/2,„^®4/2„^. (163) 

mi ,mi 

This operator acts on the direct product of two Hilbert spaces of spins at site 
1 and site L. It has the property that 

*J^^^|vac)i (g) |vac)L = | J, M)i,l. (164) 

If we define a set of degenerate VBS states {\VBSl{J, M))} such that 

\VBSl{J,M)) = ^'1^j\YBSl), J = 0,...,S, M = ~J,...,J, (165) 

then these (5* + 1)^ states (|165p are not only linearly independent but also mutually 
orthogonal. 



5.3.3. The Orthogonality 

To show the orthogonality of the degenerate VBS states (|165p . it is convenient to 
introduce the total spin operators of the subsystem (block) : 

L L L 

^t+ot-E«]^^' S^,,=J2b]a„ S^,,^-Y,{a]a,-b]b,). (166) 
j=i i=i j=i 

First we show that the set of operators {5*^^)., S^^^., S^^^} commute with the product 
of valence bonds, i.e. 
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These commutation relations (jl67p can be shown in similar ways. Take the commu- 
tator with S^^^ first. We re- write the commutator as 

L-l 

[^t+ot, n(«]^Vi-^H+i)''] (168) 

= J2i<^\bl ~ h\o}^^ . . . [ 5+ + S;+ 1 , (a]6]+i - h]a]^^f ] • • • 

Then using commutators [aj,aj] = bij and — Sij, we find that 

= [a% + , (a]bl+i - bla]^if] 

= a%A'^]b]^l-b]al,f]+al,[b,+^,ia]b]^,^b]al^,f] 

= 0. (169) 
Therefore [S^^^ , nt=i^(a]^]+i " ^l^l+i)^] = 0- In (UMl) we have used 

In a similar way, we find that the commutator with S^^Q^ also vanishes. Next we 
consider the commutator with S^^^: 

[^tot, nVjfei+i-&H+i)''] (171) 

= E(«!4 - b{alf ■■■[S^ + SUi. - &la]+i)^] • • • 

•••(ai_i6[-6Liai)^. 
In the right hand side of (|17ip . the commutator involved also vanishes because 

[^ + ^i,(a]&l+i-&]a]+i)^] 



2' 

'-alia, , {a]bl, - fejat^J^] - h% , (a^bl^, - ^ja]^,)^] 

(172) 
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Substituting into p7T|) . we obtain [S^^^ , Hj^r/CaJ^i+i - &]a]+J-^ ] = 0. Now 

we shall show that the state |VBSi,(J, M)) is a common eigenstate of S^^^ and the 
total spin square S^^^ = ^{S^otS^ot + '^tot'^tot) + ('5'tot)^ with eigenvalues M and 
J{J + 1), respectively. Using the commutation relations (|167p . we can show that 

L-l 

L-l 

S^,JVBSl{J,M)) = l[{a]b]+, - bla]^,f{S'i + J, Af)i,L|vac)2,...,L-i. 

(173) 

Then from the definition of the state |VBSl(J, M)) and the following relations: 

[St + StM M)i,L = ^iJTM){J±M + l) I J, M ± 1), 
(St + Sl)\J,M)i^L = M|J,M)i,L, (174) 

we obtain 

5ttt|VBSL(J,M)) = ^{JTM)iJ±M + l) |VBSl(J,M±1)), 
S^,,\YBSLiJ,M)) = M\YBSl{J,M)) (175) 

and hence S'?„JVBSz,( J, Af)) = J(J + 1)|VBSl(J, M)). It is now proved that 
|VBSl(J, M)) is a common eigenstate of Hermitian operators S^^^ and S'^^t with 
eigenvalues M and J(J + 1), respectively. Therefore the states with different eigen- 
values (J, M) are orthogonal to each other. Thus we have proved the orthogonality 
of the set {|VBSl(J,M)) | J = 0, . . . , S"; A/ = - J, . . . , J}. 



5.3.4. Completeness and Equivalence 

It is obvious from (jl62p that any ground state |G; J, ft) can be written as a linear 
superposition over these degenerate VBS states: 

\G;J,Cl) = \[^^±^±^l^^^^ E Xjm\YBSl{J,M)), (176) 



and vice versa. Now we can derive the completeness relation of the set {|G; J,Cl)} 
using ([TM|) . (HH]) and ([TMl) : 

df2|G; J,i7)(G; J,J1| (177) 
4. (S + .H-1,;(5-,;)^ ^ *t^,|VBS.,(VBS.|*.„^ 



(2J+1)! 2,7 + 1 



The set of states {\['jjy^|VBSL) | = —J,..., J} are linearly independent. So 
that the rank of {|G; J, (i)} with fixed J value is 2J + 1, which can be obtained 
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from the completeness relation (I177P (see ^33) . Thus the total number of linearly 
independent states of the set {|G; J, 17)} is ^j^qC^J + 1) = (5* + 1)'^, which is 
exactly the degeneracy of the ground states of (|156p . So that {|G; J, fj)} forms 
a complete set of zero-energy ground states. The set {|VBSl( J, M))} differs from 
{|G; J, f2)} by a change of basis, therefore it also forms a complete set of zero-energy 
ground states. These two sets ()158p and (I165|) are equivalent in description of the 
degenerate ground states of the block Hamiltonian (I156p . (More details such as the 
expansion ()162p etc. can be found in Chapter 9 of l^.) 

5.4. Eigenvectors of the Density Matrix 

Eigenvalues of the density matrix Or can be derived indirectly, as in^^ for spin-1 
(see ? I4.4I for comparison) and in for spin-S*. The basic idea is the following: 
Because the density matrix is independent of both the total length of the spin chain 
and the starting site of the bl ock w e can add boundary spins directly to the ends 
of the block. It was shown in by a Schmidt decomposition (see Section 2.5 

ofEHl) that non-zero eigenvalues of the density matrix (|153l) are equal to those of 
the density matrix of the two boundary spins. All other eigenvalues of the density 
matrix (jl53p are zero. This fact reveals the structure of the density matrix as a 
projector (up to a multiplicative 'scaling' matrix) onto a subspace of dimension 

Now we propose a theorem on the eigenvectors of the density matrix given by 
(|153p . The explicit construction of eigenvectors allows us to diagonalize the density 
matrix directly. The set of eigenvectors also spans the subspace that the density 
matrix projects onto. 

Theorem 3. Eigenvectors of the density matrix U53\) with non-zero eigenvalues 
are given by the set { |G; J, f}) } il58\) . or, equivalently, by the set { |VBSl(J, M)) } 
I1165\) . i.e. They are zero-energy ground states of the block Hamiltonian Hb I1156\) . 
(We shall emphasize that the eigenvectors \G;J,Q) and |VBSl(J, Af)) correspond 
to the vectors |A) in ^3.1.1\ and in Theorem[B of ^ [KM . ) 

Proof: We prove this theorem by showing that the density matrix (jl53l) can 
be written as a projector in diagonal form onto the orthogonal degenerate VBS 
states { |VBSl(J, M)) } introduced in (I165p . An alternative proof taking a different 
approach is given in the next section §[53]. 

First, it is realized from the definition of spinor coordinates ()134p that if we 
change variables (m, u) to {iv* , —iu*), then the unit vector D, is inverted about the 
origin to —fl. So that we have ■° 

{u*b^ -v*a^f\Yac) =i'^Vs^.\~n), (178) 

where | — f2) means a spin-S'/2 coherent state for a point opposite to ft on the unit 
sphere. Therefore, taking expressions of the boundary operator (|154p and the 



June 11, 2010 0:33 WSPC/INSTRUCTION FILE AKLT 



42 V. E. Korepin & Y. Xu 



block VBS State |VBSl) we have 

StlVBSi) = 

s 



(179) 



L-l 



S\ Jl (a]^]+i ~ I - ^^o>i «> |vac)2 «) • • • » |vac) 



L-l 



L + IIL- 



Consequently the density matrix (llSSp can be re-written as 



Pl 



S + l 
(25+1)! 



S\S\ 



L-l 



5* 



(180) 



L-l 



■I 



5+1 



I vac; 2 (vac I 



|vac;i_i(vac| 



4+1 n 



where and are (S* + l)-dimensional identities associated with site 1 and 
site L, respectively. In obtaining (jlSOl) . we have changed integral variables from Q,q 
, f^L+i to —CIq, — r2i_|_i and performed these two integrals using the completeness 
relation (I136p . Next we notice that (see 

s J 

4+1^ 4+1 = E E \j,m),,l{j,m\ (181) 

J=OM=-J 
S J 

= E E *jMlvac)i(vac| «) |vac)L(vac|*,7A/. 

J=OM=-J 

As a result, combining (|180p and (jlSip . recalling definitions of |VBSl) (|155p and 
|VBSl(J, M)) (|165p . the density matrix takes the following final form 



Pl = 



S + 1 



1 L 



(25+1)! 
S + l 



S\S\ 



5- 



lE E *!7Af|VBSL)(VBSi|vI' 



JM 



(182) 



(25+1)! 



5!5! 



J=0 M=-J 
J 



lE E |VBSi(J,Af))(VBSi(J,M)|. 



J=0 Af=-J 

The set of degenerate VBS states {|VBSl( J, Af))} with J = 0, . . . , 5 and M = 
— J, . . . , J forms an orthogonal basis (see §[53]). These (5 + 1)^ states also forms 
a complete set of zero-energy ground states of the block Hamiltonian (|156p (see 
? l5.3n . So that in expression (|182l) we have put the density matrix as a projector in 
diagonal form over an orthogonal basis. Each degenerate VBS state |VBSl(J, M)) 
is an eigenvector of the density matrix, so as any of the state |G; J, fi) (because 
of the degeneracy of corresponding eigenvalues of the density matrix, see ^i lS.GI and 
§ 15.71 that the eigenvalues depend only on J). Thus we have proved Theorem [3l □ 



5.5. An Alternative Proof of Theorem\^ 



It was shown in i) l5.4l that the density matrix takes a diagonal form in the basis of 
zero-energy ground states of the block Hamiltonian (|156p . In this section, we show 
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the same result by taking a different approach. This alternative proof of Theorem 
[3] does not involve the coherent state basis. 

The proof uses the fact that the density matrix is independent of the starting 
site and the total length of the chain (see ^^ 15.21 ). So that we could change the 
configuration of the whole system by adding the two ending spins directly to the 
block without affecting the form of the block density matrix. The new system now 
has L + 2 sites with the block starting at site 1 and ending at site L. Let us start 
with the ground VBS state of the Hamiltonian ()130|) of the system with N ^ L: 

L g 

|VBS) ^ n - l^ac)- (183) 

In order to calculate the density matrix — tro,L+i [p], where p is defined in 
(|133p . we introduce a useful identity: 

(_l\S-J+M 

0,L + l(J,M|(|s)o,l®k)L,L + l) = ^ ^^^^^ |J,-M)i,L, (184) 

where |J, M)o_l+i is identical to the spin state defined in (|164p except for site 
indices. \s)ij in (I184p is the normalized singlet state with S valence bonds defined 
as 



(-1) 



S 5/2 



= J2 (-in5/2,-m),® |5/2,m)j. (185) 



^ • - m=-SI1 

Identity (|184p is derived using properties of the singlet state (|185p and Clebsch- 
Gordan cocfhcients as follows: 

o,L+i(J,M| (|s)o,i ® 

^ (J,M|^/2,mo;5/2, 

.i-=4L (-irM^/2,-mi)o|^/2,mi)i 

V'-' + t , 

mi,= — S/2 



E (-l)'"°+'"^+U^,M|^/2,mo;5/2,mi+i) 

■\S/2,-mo)^\Sl2,-mL+i)L- (186) 
Here the Clebsch-Gordan coefficient is defined by 

(J,M|^/2,TOo;^/2,TOL+i) = (|^/2,TOo)« ® \S/2,mL+i)j) . (187) 
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Then using the symmetry property of Clebsch-Gordan coefficients 

(J,M|5/2,mo;5/2,mi+i) = i-lf-''iJ,^M\S/2,-mo;S/2,-mL+i) 

(188) 

and the completeness of the basis { \S/2,mo)o \S/2,mL+i)L+i }, we obtain the 
identity (jTMl) . 

With the help of identity (|184p . we calculate the partial inner product of the VBS 
state (|183p with the state | J, A/)o,l+i, which is involved in taking trace of boundary 
spins. The VBS state |VBS) is decomposed into the bulk part and boundary parts, 
then making use of (|184p . we have 

o,L+i(J,M|VBS> 
= o,L+i(J, M\ Yl (a]&]+i - 6]a]+i) |vac) 

3=0 

L-1 g 
= Sl{S+iy. n - ^Wj + i) 0,L+l(J,M|s)o,lk>L,L + l|vac)2...L-l 

= (Slf fl - b]al,y i-lf-'+'''\J, -Af)i,,.|vac)2...i_i 

= {-lf--'+^\S\)^\YBSLiJ,~M)). (189) 

We see that the {S + 1)^ degenerate VBS states { |VBSl( J, M)) } defined in (ITC5|) 
appear in the partial inner product (jl89p . As discussed in S I5.3I . they form a complete 
set of zero-energy ground states of the block Hamiltonian (|156l) . 
Now, it is straightforward to evaluate the density matrix as 

o,L+i( J, M|VBS) (VBS| J, M)o,L+i 



tro,L+i [p] = 



(VBS|VBS) 

^^^^||^5:|VBS.(J,-M))(VBS.(J,-M)|. 



(51^4 



(190) 

This expression is identical to (jl82p as we change dummy index from M to —M. 
Therefore, in this approach again we arrive at Theorem[3]that the density matrix is 
proportional to a projector onto a subspace spanned by the {S + 1)^ ground states 
of the block Hamiltonian (|156|) . Normalization (VBSjVBS) has been obtained in 
§ 15.1.31 . States |VBSl(J, M)) have been shown to be mutually orthogonal in 



5.6. Eigenvalues of the Density Matrix (normalization of 
degenerate VBS states) 

As the next step in analyzing the spectrum of the density matrix, now we study the 
eigenvalues. Based on the diagonalized form (jl82|) . it is clear that eigenvalues of the 
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density matrix can be derived from the normalization of degenerate VBS states. 
We obtain an explicit expression for eigenvalues in terms of Wigner 3j-symbols in 
this section. 

First, the following property is important: Normalization of the degenerate VBS 
state |VBSl(J, M)) depends only on J and is independent of M. This point is 
important in proving that any |G; J, 17) is an eigenvector of because it can be 
written as a superposition of |VBSl( J, Af ))'s with the same J value (|176p . With 
the introduction of total spin operators of the block S^^, S^^^ and Sl, (see §0), 
we prove the statement as follows: 



(VBSl(J, Af ± 1)|VBSl(J,M± 1)) 
1 



(J T M){J ± M + 1) 
1 

(J T M){J ± M + 1) 
(VBSl(J,A//)|VBSl(J,A/)). 



(VBSi( J, M)\S^,,St,\VBSL{J, M)) 



(VBSi( J, M)\{St,, - {Sl,Y T 5t^ot)|VBSi(J, M)) 



(191) 



Here we have used the fact that |VBSl(J, Af)) is the common eigenstate of S^^ 



and with eigenvalues J( J + 1) and Af, respectively (see §[0]). 

It is also realized that the normalization of |VBSl( J, M)) can be calculated from 
integrating the inner product of |G; J, Cl) with itself over the unit vector (l such that 



1 

Ait 



dJ7(G; J,f2|G; J,f2) 



{S + J + 1)\{S - J)\J\J\ 
(2./+1)! 



(VBSl(J,A/)|VBSl(J,A/)). 



(192) 



In obtaining this relation (|192l) we have used expansion p76p and orthogonality 
(fTCDD in ^ [Ql . 

Let us consider the integral involved in (|192p . Using the coherent state basis 
(|135p and completeness relation (|136|) as before, we obtain 



in ' 



(193) 



1 

An 



(25+1)! 



1 L 



47r 



S-J 



L 



L-l 

n 



-(i+rii-fi) 



-(1 - ■ flj+i] 



-{1 + n-nL) 



Now we expand ^{1 ~ Cli ■ Clj) in terms of spherical harmonics as in (|145p . then 



June 11, 2010 0:33 WSPC/INSTRUCTION FILE AKLT 



46 V. E. Korepin & Y. Xu 

integrate over tl and from to ^Il-i, the right hand side of (|193|) is equal to 
44(25 + 1)!)^ VVV V V V 

{s + 1)1^-^3 - J + + ^ ^ ^ 

dQi J dClL X^^\li, S)X{Il, S - J)X^{1, J) 
Here we apply the following useful formula: 



{2h + l){2lL + l)i2l + l) fhlL l\ f h II I 



An \ 0/ \mi niL m 

(195) 

where [ ^ ^ I is the Wigner 3j-symbol. Using formula (|195p . we carry out 
\mi rriL m J 

the remaining integrals in (I194p and obtain 



(5 + l)i-i(5-J + l)(J + l)2 



((2^+1)!)^ 



EEE E E E 

li—Olj^—O l—O mi— — /i mL — — lL m— — l 



{2h + 1){21l + 1K21 + l)X''-\h,S)X{lL,S- J)X'{1, J) 



\ 2 / ^ 2 



y \ mi rriL m ^ 
These 3j-symbols obey the following orthogonality relation: 



(196) 



Y: i2l + l)(^' ^)(^' ^\)=Su,S„.rn'. (197) 

^ — ' \ mi rriL m J \ mi ra^ m I 



nil ,niL 



Using this orthogonality (|197p . we can recast expression (jl96p as 



((2g+l)!)^ ^ ^"-^ 

(5+1)^-1(5- J+1)(J+ 1)2 



E E E (198) 



{2h + l){2lL + l){2l + l)X'^~\h,S)X{lL,S-J)X^{l,J) ' 







The explicit value of ( ^ ^ 1 is given by 







h II I 




(-1) 



9 . 



'{2g-2h)\{2g-2lL)\{2g-2l)\ g\ 



(2.g+l)! {^g^l,)l[g-l^)l[g-l)V 



(199) 
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ii h -\- II + I ~ "^g (g G N), otherwise zero. Finally, the normalization of degenerate 
VBS states |VBSl(J, A/)) is obtained as 

(VBSl(J,M)|VBSl(J,M)) (200) 

(2J+1)!((25+1)!)^ ^ ^ 



{S+l)^-HS + J+l)l{S-J+mJ+l)l{J+l)\ 



X! X! X! 



{2h + + + l)A^-i(?i, S)\{Il, S - J)A^(?, J) 



h II I 




Combining results of (|182l) and (|200p . we arrive at the following theorem on 
eigenvalues: 

Theorem 4. Eigenvalues A( J) ( J = 0, . . . , 5) o/ the density matrix pj^ are inde- 
pendent of D, and M in defining eigenvectors (see il58]) and il65]) ). An explicit 
expression is given by the following triple sum 

A(J) (201) 



S + 1 



1 L 



(25+1)! 



S\S\ 

^(VBSi( J, M)\YBSl{J, M)) 



(2J + 1)!5'!6'! ^ t-.^ ^ 

- (5 + J + ms - J + 1)!(J + 1)!(J + 1)! ^^^^h 

{2h + 1){21l + 1)(2/ + l)A^-i(?i, S)\{Il, S - J)X'{1, J) (^^J J 

Although not straightforward to verify, this expression (|20ip should be consistent 
with eigenvalues given through the recurrence expression (j214p in the next section 
it \5.7\ and the expression Aq, in ii l4.4l as a special case. We could check the case when 
5 = 1 that 

(VBSl(0,0)|VBSl(0,0)) ^ ^(3^ + 3(-l)^), 

(VBSl(1, A/)|VBSi(l, M)) = i(3^ - (-1)^), (202) 

where we have used the selection rule of the Wigner 3j-symbol. From ()148p we find 
that (VBSjVBS) — 2 • 3^, so that we obtain the correct eigenvalues of the density 
matrix from the above result pOip (see § 14.41 for comparison). 

We shall emphasize at this point that given eigenvalues pOip , both von Neumann 
entropy 

s 

N = -tr [ Pi In ] - ^ (2 J + 1) A( J) In A( J) (203) 

j=o 

and Renyi entropy 



Sr. = ^ In {tr [ p2 ] } = ^ In I ^ (2 J + 1) A" (J) | 



(204) 
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can be expressed directly. 

5.7. Eigenvalues of the Density Matrix (recurrence formula) 

Having constructed eigenvectors, there are more than one way to specify the cor- 
responding eigenvalues. An explicit expression of eigenvalues has been obtained in 
^ 15.61 . In this section we express eigenvalues through a conjectured recurrence for- 
mula as appeared in '25'48|_ j^g^ apply the density matrix ()153p to the state 
|G; J, (l) and obtain 



pL\G\J,n) 
S + l 



(25+1)! 



(477) 



dl7odi7L+iBt|VBSL)(VBSL|BA;f^|VBSL). 



(205) 



Using the coherent state basis (|135l) and completeness relation (|136p . the factor 
(VBSl|BA|^|VBSl) in ([205]) can be re-written as 



(VBSlIBAJ^IVBSl) 



(206) 



(25+1)! 



47r 



I (mi 



The factor ^(1 — ■ %+i) under the integral of (|206p can be expanded in terms 
of Legendre polyn omials and further in terms of spherical harmonics as discussed 
m §EX3] (see alsol^^. Using the expansion (I145P and orthogonality of spherical 
harmonics, the integrals over Ctj with j = 2, . . . , L — 1 in (|206p can be performed. 
The result is 



(VBSl\BAI\YBSl} 



S + l 
(47r)2 



(25+1)! 



L S 



1=0 



(WlUi-UlMi) (WW^+W^) (UlVl+1 - vlul+i) ■ 



(207) 



We plug the expression ()207p into ([205]). Using transformation properties under 
SU{2) and a binomial expansion (see S I5.8n . the integral over (Iq yields that 



da 



(u'^b\ - v^a\^ {uqVi - vqUi) 



An 



J (^uia\ + vib\ 



(208) 
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Similarly we can perform the integral over f^L+i- As a result, the following expres- 
sion is obtained from (|205p : 

1 ^ 

pJG; J, n) = — ^ ^(2/ + l)A^-i(;, S)Kl{n) IVBSl) . (209) 

^ ' 1=0 

The operator Kj{tl) involved in (I209|) is defined as 

■ (uul + vvlY (^ULa\ + VLb^j}j ^ Pi{hi ■ Ql). (210) 

It is expressed as an integral depending on the order I of the Legendre polynomial 
Pi{Cli - CIl). Kj{tl) can be calculated from the lowest few orders (see ^ \5.8\ for exam- 
ple). It becomes increasingly difficult to perform the integral as order I inc reases. 
Based on the eigenvalue expressions of the density matrix obtained inmSi, we 
make a conjecture on the explicit form of the operator Kj{tl) for generic order /: 
Conjecture 1 

Kjm - (^) ' // J( J + 1) - Isi^S + 1)) 4. (211) 

Here the polynomial // (x) satisfy the recurrence relation 

21 + 1 f Ax A . . . I fS-l + lV ^ , , 

j Ii-i{x) 

(212) 



with Iq = 1 and Ii = j^frj^- 

Note that it is important that Kj{Ct) (x Aj defined in ()157p and Ii{x) has the 
same order as the Legendre polynomial Pi{x). The recurrence relation pi2p was 



proposed in and used in HSl to obtain the eigenvalues of the density matrix. 
(The original definition of Ii{x) differs from our definition in (|212p by a factor of 
{21 -\- 1)/Att.) Conjecture 1 is an alternative form of Theorem [3] together with 
Theorem 21 which also gives eigenvalues through the recurrence relation (|212p . 
Indeed, expressions (|209l) . altogether with (I21ip and (|212l) yields that 

PL\G;J,n) (213) 
- ^ ^{2l + l)X'^-\l,S)Ii (Ij{J + 1)-1:S{Is+1)) \G;J,n). 



(^+1)2 ^' ^ ' V2 ' '2 ^2 



Non-zero eigenvalues ( J = 0, 1, . . . , S*) are seen from (I213P as 



(214) 
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Since all other eigenvalues of the density matrix are vanishing, then we conclude 
again that the density matrix (I153P is a projector onto a subspace of dimension 
{S + 1)^. This subspace is spanned by the set of vectors {|G; J, fi)} (I158p . (The 
rank of the set is equal to (S' + l)^.) Furthermore, we observe from (|214p again that 
non-zero eigenvalues A(J) depend only on J, not on tl. Therefore, {|G; J, O)} with 
fixed J value spans a degenerate subspace with the same eigenvalue. 



5.8. The Large Block Limit 

In th e limit L — )■ oo, that is when the size of the block becomes large, we learned 
fi-Qjn | 18|33|48 | von Neumann entropy reaches the saturated value S'v n = 

In (S" + 1) . This fact implies that the density matrix (denoted by in the limit) 
can only take the form (see for a general proof) 



$oo, (215) 



where /(5_|_i)2 is the identity of dimension (S'+l)^ and $oo is an infinite dimensional 
matrix with only zero entries. In this section, we give a proof of Conjecture 1 (|21ip 
in the limiting case as L — > oo. Then we shall verify the structure of the density 
matrix (I215|) explicitly. 

We first realize from p46p that as L — > oo, X^^^{l,S) — > (5;_o- Therefore only 
the first term with Z = is left in ()209p and contributes to the final result. So that 
we need only to calculate Kq{Ci): 

s 



■ {ulvl - vlulf-' {uul + vvlf (uLal + VLblY . (216) 

It is useful to know transformation properties of the integrand in (|216p under 
SU{2). The pair of variables {u, v) defined in (|134p and bosonic annihilation opera- 
tors (a, 5) in the Schwinger representation both transform as spinors under SU{2). 
That is to say, if we take an arbitrary element D e SU{2) (a 2 x 2 unitary matrix 
with determinant 1), then {u,v), etc. transform according to 

On the other hand, {u*,v*), {~v,u), {a^,b^) and (—6, a) transform conjugately to 
{u,v). That is to say {u*,v*), etc. transform according to 

:•>-(;')■ (-) 

The combinations appeared in K^{tl) (|216p 

uia\+vib\, uul+vvl, ulvl — vlu*j^, uu*j^+vv1, ULa]^ + VLb\ (219) 
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as well as in (|157p . boundary operator in (|154p . etc. all transform covari- 
antly under SU{2), i.e. those expressions keep their form in the new (transformed) 
coordinates. 

These transformation properties (|217p . (|218l) can be used to simplify the Kq{Ci) 
integral. We first make a SU{2) transform 

under the part of the integral ()216p over Cti. Then this part of integral becomes 

dCli (^uia\-hvib\y (uul + vviy^ {^vlf'\ (221) 

This can be calculated using binomial expansion and the result is 

J^(^ua\+vb\y [~b\y~' . (222) 

Then we make an inverse transform in (j222l) using D^^^ = ' consequently (I216P 
is put in a form with a single integral over Cl^ remaining: 

47r / + .i\J 



Kl{n)^j^{ua\+vbl) (223) 
■ J dClL (^a\vl ~ b\ulj {uul-\-vviy-'(uLa^^-\-VLb'^j^^ . 
Now we make another SU(2) transform using 

^u(:) = (i), (224) 



-V u J V^/ \U, 

then the remaining integral over in (|223p becomes 

dQL [a\vl - b\ul^ ^"^ (ulf (uLa[ + ^ . (225) 

Using again binomial expansion, this integral (j225l) yields 

4^ ^t^,t K^jV'^fj^-^ 



^i^albl-blal) (^a[j . (226) 



At last we make an inverse transform in (j226p using D„ ^ and plug the result 

into (|223p . the final form is 

Klm^l^-^y A\. (227) 

This expression is consistent with Conjecture 1 (|21ip . which also proves that 
{|G; J, is a set of eigenvectors of the density matrix as L — >■ cx). Let us denote 
the density matrix in the limit by p^. Then ([227| leads to the result (see (|2T3l) ) 

p^\G;J,n) = ^^^\G;J,n). (228) 
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We find from l|228p that the hmiting eigenvalue 

Aoo - (^^Yj2 ' ^ ^ °° (229) 

is independent of J. Any vector of the {S + 1) ^-dimensional subspace spanned by 
the set {\G;J,fl)} is an eigenvector of with the same eigenvalue 1/(5' + 1)^. 
Therefore acts on this subspace as (proportional to) the identity /(5+i)2. So 
that we have proved explicitly that the density matrix takes the form (|215p in the 
large block limit. The limiting density matrix is proportional to a projector P(5_|_x)2 
on the degenerate ground states of the block Hamiltonian Hb (|156p 

Pl ^ Poo = (^sll)'^ P(s+i)= ' L^oo. (230) 

In addition, we also derive from the eigenvalues that the von Neumann entropy 
'S'v N = ~X]j=o(2'^+ 1) Aoo In Aoo coincides with the Renyi entropy S'Fi(a) = 
In 7=o(2^ + 1)A^| and is equal to the saturated value 

N = 5'R,(a) = ln(5 + l)^ L ^ oo. (231) 



5.9. The Density Matrix and Correlation Functions 

The relation between the density matrix and correlation functions was studied in 
[5 , 46 | 48 i 76] ^ .^g^g gjjown in that the density matrix contains information of all 
correlation functions in the ground state. The original proof was for spin S — 1/2. 
In this section we generalize the result to generic spin-S" as in and the proof is 
written in a form applicable but not restricted to the VBS state. 

The Hilbert space associated with a spin-S' is (25' + l)-dimensional. Therefore 
we could choose a basis of (25 + 1)^ linearly independent matrices such that an 
arbitrary operator defined in the Hilbert space can be written as a superposition 
over the basis. Let us denote the basis by {Aat | a, 6 = 1, . . . , 25"+ 1}, in which each 
matrix Aat is labeled by a pair of indices a and b with totally (25" + 1)^ possible 
combinations. The matrix element is defined as 

iAab)ki ^ 6ak6bi, kJ = l,...,2S + l. (232) 

In addition to {Aaf,}, we introduce an equivalent conjugate basis {Aab} such that 

{Aat)ki ^ SaiSbk, a,6,fc,/ = 1,...,25 + 1. (233) 

These matrices (|232p and ()233p are actually matrix representation of operators 
{\S, to) (5, m'l I m,m' — —S, . . . , S}. They are normalized such that 

tT{AabAcd) = ^^{Aab)kl{Acd)lk = ^ SalSbk^d^dk = SacSbd- (234) 
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Here 'tr' takes trace at one and the same site. Because of the completeness of {Aab} 
at each site, the density matrix of the block (B) can be written as (see (|96|) ) 



Ps = tTE\G){G\ = I (^^a,fc, I coefr{a,5,}, 

{ajbj} \jeB 



(235) 



where |G) denotes the unique ground state {e.g. the VBS state), trg takes traces 
of sites outside the block (i.e. sites in the environment E) and coeS{ajbj} denotes 
the coefficient. Using the normalization property (I234p . the coefficient coeS{ajbj} 
with label j taking values within the block can be expressed as 

coeS{ajbj} ^ ^ J| tr(Aa^b^.v4c,<i,)coeff{cjdj} 
{cjdj}jeB 




Aa,b, |G)(G| 



(G| I (^Aa,i>, I |G). 



(236) 



Here tr^ takes traces of sites within the block and tr takes traces of all lattice sites. 
Combing p35p with (|236p . we have the final form 



Pb 




(237) 



This is the expression of the density matrix with entries related to multi-point 
correlation functions (G| {^j^g Aa^b,^ |G) in the ground state. All possible combi- 
nations {cLjhj} are involved in the summation. Therefore, we have prove for generic 
spin-S* that the density matrix contains information of all correlation functions. The 
matrix elements are all multi-point correlators. 



6. The One dimensional Inhomogeneous Model 

The most general 1-dimensional model is the inhomogeneous model in which spins 
at different lattice site can take different values. As a special case of the generalized 
model defined in ii l2.4l . we associate a positive integer number (called multiplicity 
numbers, see ii l2.4.1l ) to each bond of the lattice and denote by (M^- = Mji) 
the multiplicity number between sites i and j. They are related to bulk spins by 
the following relation which ensures the existence of a unique ground state 



(238) 
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with 25*0 = A/qi and 2Sm+i = Mat^at+i for ending spins. (Equation (j238p is a special 
case of the more general relation ([30]).) The condition for solvability of relation (|238p 
is 

N+l 

Y^{-iyS,=Q. (239) 
Solution to relation (|238p under condition (|239p is 

M,-,+i = 2^(-l)^-'5,>l. (240) 

1=0 

(More details can be found inl^.) Now we defined the Hamiltonian of the inhomo- 
geneous AKLT model according to ((28)) as 

^ = E E a/(j, J- + l)^j(j,J + l). (241) 

j=0 J^Sj+Sj + i-A/j.j + i + l 

Here the projector 7rj(j,j + 1) describes interactions between neighboring spins j 
and j + 1, which projects the bond spin Jjj+i = Sj + Sj+i onto the subspace with 
total spin J {J = Sj + Sj+i — Af,-,,+i + 1, . . . , Sj + Sj+i). An explicit expression of 
TTj{j,j + l) is given in ^ \2.2\ and The coefficient Cj{j,j + 1) can take an arbitrary 
positive value. This Hamiltonian p4ip has a unique ground state (VBS state, see 

§■231) • 

Following L^-iJ, we study the entanglement of the unique VBS ground state of the 
inhomogeneous model (j24ip in this section. 

6.1. The VBS Ground State 

The unique VBS ground state of the Hamiltonian (|24ip is given in the Schwinger 
representation by 

N ^ 

I VBS) ^ n (a]6]+i - 6]a]+i) |vac), (242) 

where a\ W are bosonic creation operators defined in exactly the same way as in 
? I2.4.2I . the constraint on the total boson occupation number is now (ajaj+6]&j)/2 — 
Sj. The pure state density matrix of the VBS ground state (|242p is 

|VBS)(VBS| 
P= (VBSIVBS)- 

Normalization of the VBS state is (calculation similar to those in ? l5.1.3n 

N+l 

n(25, + l)! 

(VBSIVBS) = -j^ . (244) 

n(M,-,+i + 1) 

j=0 
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(See 1^ for more details.) 



6.2. The Block Density Matrix 



We take a block of L contiguous bulk spins as a subsystem, which starts from site 
k and ends at site k + L — 1. Using the coherent state basis (|135p and completeness 
relation (|136p . tracing out degrees of freedom outside the block, can be written 



as^ 



77 



Pl = 



k-l N+1 

n n d% 

j=0 j=k+L 



k-2 N 

n n 

j=0 j=k+L 



-(1 - ^Ij ■ VLj + i] 



Mi 



(245) 

B^|VBSl)(VBSl|B 



fe+L-l 

n 

j=k 



i^S, + 1)! 
47r 





Af+l 


N 


7 




n 




j=0 


J=0 



Here the boundary operator B and block VBS state |VBSl) are defined as 
B = {uk-ibk - Vk-iak)^''"-^'" {ttk+L^iVk+L - bk+L-iUk+h)' 



k+L-2 

YBS,)^ n {-]bli-b]ali) |vac), 



(246) 
(247) 



respectively. After performing integrals over ttj {j = 0, 1, . . . , fc — 2, fc + L + 
1, . . . , A'', A'^ + 1) in the numerator and all integrals in the denominator, the density 
matrix turns out to be independent of spins outside the block. This Prop erty has 
been proved for the homogeneous AKLT model in §[0] (see also Ll^l^^Lmi) . There- 
fore, we can re-label spins within the block for notational convenience. Let fc = 1 
and the density matrix takes the form 

L 

J dnodnL+iB^VBSL){yBSL\B (248) 

n(2^, + l)! 



Pl = 



with 



Mr 



L-1 



Mi, 



I vac). 



(249) 
(250) 



Again, the remaining two integrals in (12481) are kept in the present form for later 
use. 
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6.3. Ground States of the Block Hamiltonian 

The block Hamiltonian with the re-labeling fc = 1 reads 

L — 1 Sj+Sj + \ 

Hb = Y. E Cj{j,j + l)Pj{j,j + l). (251) 

j=l J=Sj+Sj + i-Mjj + i+l 

Now the degenerate ground states are constructed in a similar way as in S I5.3I . The 
new A^j operator is defined as: 

A\ EE (ua\ + vb\) (al6[ - b\al) {ual + , (252) 

where J_ = (Mqi - Ml.l+i)/2, J+ ee {Mqi + Ml,l+i)/2 and \ J-\ < J < J+. Then 
the set of ground states of the block Hamiltonian (|251l) is 

\G;J,Ci)=AI\YBSl), J =\J-l...,J+. (253) 

To prove that any state |G; J, f2) is a zero-energy ground state of (|25ip . we essentially 
repeat the arguments as in ^ \b.'3\ for any site j and bond + 1): 

(1) The total power of and 5| is 25^, so that we have spin-5j at site j; 

(2) -i(Af,-i,, + Af,+i,,+2) < Jl,+, ^ ^1 + ^1+1 < + Ai,+i,,+2) by 
a binomial expansion, so that the maximum value of the bond spin Jjj+i is 
^(Mj_ij -I- Mj_|_ij_|_2) = Sj + Sj+i — Mjj+i (from SU{2) invariance, seel^. 

Therefore, the state |G; J, Cl) defined in (|253p has spin-S'j at site j and no projection 
onto the Jjj+i > Sj + Sj+i — Aljj+i subspace for any bond. 

Parallelly, we also introduce an orthogonal basis in description of the degenerate 
zero-energy ground states of Hb (|25ip . i.e. the degenerate VBS states. Using the 
sam e nota tions as in § 15.31 . the operator defined in p52p can be expanded as 
(seeElIIZl) 

,t / (J+ + J + 1)!(J- + J)!(,/+ - J)!(- J- + J)! ^ 

= y ^j—^ XjM (254) 

mi+m_L=_A/ ^ ^ 
mi ,mL 

where (iMoi, toi; ^Ml.l+i, 'mL\J: M) is the Clebsch-Gordan coefficient. Again, the 
particular form of the sum over mi and tol in (I254p is identified as a single spin 
state creation operator 

(255) 

So that the set of degenerate VBS states {|VBSl( J, A/))} is defined as 

|VBSl(J,A/)) = *J,^,,|VBSl), J= |J_|,..., J+, M^-J,...,J. (256) 
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Then these (Mqi + 1){Ml,l+i + 1) states (|256p are mutually orthogonal (the proof 
is exactly the same as in §[5]3]). Also, the state |G; J, ft) written as a linear super- 
position over these degenerate VBS states reads 



\G;J,{l) = 



(J+ + J + 1)!(J- + J)!( J+ - Jy.j-J- + J) 
2J+ 1 

,7 

• X7m|VBSl(J,A/)). (257) 

Therefore, as seen from (|257p . the rank of set of states {|G; J, 17)} with the 
same J value is 2 J + 1 and the total number of linearly independent states of 
the set {|G; J,f2)} is T,jt\j_\{'^J+^) = (^^oi + 1)(Ml,l+i + 1), which is exactly 
the degeneracy of the ground states of (|25ip . So that {\G;J,Cl)} forms a complete 
set of zero-energy ground states. The orthogonal set {|VBSl( J, M))} also forms a 
complete set of zero-energy ground states, which differs from {|G; J, f^)} by a change 
of basis. 

6.4. Diagonalization of the Density Matrix 



The density matrix is diagonalized in ^ \5A\ and ^ \5.5\ for the homogeneous AKLT 
model. The analysis can be made in parallel for the inhomogeneous model. 

The statement of Theorem [3] is still valid here. i.e. Eigenvectors of the density 
matrix pj^ (|248[) with non-zero eigenvalues are given by the set {|G; J, r2)} (|253p 
or {|VBSl( J, Af ))} (|256p . This exphcit construction of eigenvectors yields a direct 
diagonalization of the density matrix. 

Again, we prove the theorem by re- writing the density matrix (j248p as a pro- 
jector in diagonal form onto the orthogonal degenerate VBS states {|VBSl( J, M))} 
introduced in (|256p . 

Take expression ([248]) and integrate over Qq and using binomial expansions 
and 

da;(H-a;)™(l-x)" = ., 2"+"+^ (258) 



(rn + n + 1)! 



Then we have 



L-l 



Pl 



\ P J \ 1 



Y[i2s, + iy. 

(6l)f(4)*'^«i-f(a[)'?(6i)^^'^+i-9|VBSi) 

{YBSL\ibL)^''-'^+'-''{aLnai)^'°'-''{h)P. (259) 
The particular combinations of bosonic operators appeared in (|259l) are recognized 



up to a constant as spin creation operators ip\ . , i . , _ and 
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■i/'t ,^ 1 at site 1 and site L, respectively. Thev commute with all bond 

/ t t t t \ 

operators ( aj6j_|_j — ^joj+i ) , so that we could simplify the right hand side of 
(|259p using definition (12551) and the following identity: 

Moi + i 

E E 4M.n.iMo,-p®^lA./... + „,-iAf.,. + Jvac)i,i 

l,L(vac|ViMoi,iA./oi-p ® ^^Mr.,i. + uq-hML.L + i 
Mai ^ 

p=0 13=0 

J+ J 
= E E \J,M),,l{J,M\ 

J=\J-\M=-J 
J+ J 

= E E *!7A/|vac)i,L(vac|vl/jM. (260) 

J=|.7_| M=-J 

The resultant final form of density matrix is then 
n(2^. + l)! ./H./-|M^-. 

L-l 

= ^ Moi!Mi,i+i! ^ ^ |VBSi(J,M))(VBSi(J,M)|. 

11(25, + 1)! J=\J-\M=~j 

(261) 

The set of degenerate VBS states {|VBSl(J, M))} with J = | J_|, . . . , J+ and 
M = — J, . . . , J forms an orthogonal basis. These (A^oi + l)(Mi^L+i + 1) states also 
forms a complete set of zero-energy ground states of the block Hamiltonian pSip . 
So that in expression (j26ip we have re-written the density matrix as a projector in 
diagonal form over an orthogonal basis. Each degenerate VBS state |VBSl(J, M)) 
is an eigenvector of the density matrix, so as any of the state |G; J, Cl) (because of 
the degeneracy of corresponding eigenvalues of the density matrix, see §[6?5]) ■ Thus 
we have generalized Theorem [3] to the inhomogeneous case. 
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6.5. Eigenvalues of the Density Matrix 

Given the diagonalized form (|26ip . again eigenvalues of the density matrix pj^ are 
derived from normalization of degenerate VBS states with an explicit expression in 
terms of Wigner 3j-symbols. 

Similarly, we first calculate the integral of the norm square of |G; J, f2)over the 
unit vector fl 



— / dn{G;J,n\G;J,n) 
47r J 

{J+ + J + 1)!(J- + J)!( J+ - J)!(- J- + J)! 
(2J+1)! 



(262) 



(VBSl(J,M)|VBSl(J,M)). 



This expression (|262p also states that normalization of the degenerate VBS state is 
independent of Cl and/or M. 

Let us consider the integral involved in (|262p . Using coherent state basis (|135l) 
and completeness relation ()136p as before, we obtain 



— / dn(G;J,n\G:J,n) 
in ' 



(263) 



1 

An 



n 



(25"^ + 1) 
47r 





L 








_j=i 



J+-J 



L~l r 

n 

J-+.J 



-{I - ■ ^ij+i] 



Mi 



-J-+J 



Now we expand ^{1 — Cti ■ Ctj) 



in terms of spherical harmonics 



-1 M- ■ 1^1: 



-Y.\{i,M,,) J2 Yun,)Y;,^in,) (264) 



with 



A(Z,M,,) 



(-1)^M,,-!(M,,- + 1)! 
(M,j-0!(My +Z + 1)!' 



(265) 



Then integrate over fl and from fl2 to the right hand side of (|263p is equal 
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to 



47rn(2^.- + l)! 



L-l 



( j_ + J + 1)( j+ - J + 1)(- J- + J + 1) 



L-l 



EEEE E E 

/ /q— m— — / — — m^— — 



•A(;„, J+ - J)A(Z0, J_ + J)Xilf3, -J- + J) j dQi j dflL 

(266) 

Where we have M< = min{Mjj+i, j = 1, . . . , i — 1} and J< = min{,/_ + J, — J_ + 
J}, both being the mimmum of the correspondmg set. Now we carry out remaimng 
integrals in ([266]) using 

dCinmmYi^,mA^)Yi^^^^in) 



47r \000/\m ma mp 



(267) 



The result after integration can be further simplified by applying the following 
orthogonality relation 

/ la Ip 



E (2^/3 + 1) 



m ma rnp J \m ma mjg / p /3 p 3 



(268) 



where 



/3 



, etc. are the Wigner 3j-synibols. 

m ma mp ' 
So that finally expression (j266l) is equal to 

L 

n(25, + l)! 



L-l 



Af< J+-J ,/< 

EE E 



•(2; + i)(2;„ + i)(2Z/3 + i) 



(j_ + J + 1)( j+ - J + 1)(- J- + J + 1) 

A(/c., J+- J)A(;^, J_ + J)A(/^,-J_ + J) 

(269) 



L-l 

n A(;,M„-+i) 



ODD 
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The explicit value q ) given by 



I la Ip 




(270) 



(-1) 



g l {2g-2l)\{2g-2la)\{2g-2l0)\ 



(2ff+l)! 



[g-l)\[g-la)\{g-lp)V 



ii I -\- la + 1/3 = 2g {g £ N), otherwise zero. 

Combining results of (|261l) . (|262p and (|269l) . we arrive at the following result 
for eigenvalues: Eigenvalues A(J) (J = | J_|, . . . , J+) of the density matrix are inde- 
pendent of and/or M in defining eigenvectors (see (|253|) and (|256|) ). An explicit 
expression is given by the following triple sum 



L-l 



n(Mw+i+i) 



A(J) 



-MoiIMl,l+iI (VBSl(J, Af)|VBSL(J, M)) 



n(2^, + l)! 

(2J+l)!Moi!Mi,L+i! 



(J+ + J + 1)!( J- + J + 1)!( J+ - J + 1)!(- J- + J + 1)! 



Af< J+-./ ,/< 

EE E 



L-1 



A(;„, J+ - J)\{lp,J- + J)X{lp,-J^ + J) 



•(2; + i)(2/„ + i)(2;^ + i) 







(271) 



6.6. The Large Block Limit 

In this section, we generalize the characteristic properties (ii l5.8l ) of the limiting 
density matrix to the inhomogeneous model. 

Let us apply the density matrix pj^ (|248p to the state |G; J, Cl) (|253p and get 

PL\G:J,h) (272) 

L 

n(Af„+i+i) 

= J df2odf2i+iBt|VBSL)(VBSL|B4|VBSL). 

n(25, + l)! 

3=1 



Using the coherent state basis ()135p and completeness relation (|136p . the factor 
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(VBSl|BA|^|VBSl) in (PUSI) can be re- written as 
(VBSl|SA^,|VBSl) 



(273) 



(25-, + 1)! 



47r 



L \ L-1 

n^^. n 



-(1 - VLj ■ 



n 



Mi 



J+-J 



We plug the expression (I273P into (|272p . Using transformation properties under 
SU (2) and a binomial expansion, the integral over VIq yields that 



47r 



01 



(274) 



Similarly we can perform the integral over f^L+i. Then using expansion (|264p and 
orthogonality of spherical harmonics, other integrals over f2j with j = 2, . . . , L — 1 
in p73p can be performed. As a result, the following expression is obtained from 



(47r)2 ^ 



(2/ + 1) 



L-1 



i^/(l})|VBSL) 



(275) 



The operator k\{^) involved in (|275l) is defined as 



Mo 



J-+J 
A/i.i+i 



(276) 



It is expressed as an integral depending on the order / of the Legendre polynomial 

Plitli-ClL). 

There was no ambiguity in defining the large block limit in the homogeneous 
AKLT model (see ii l5.8l ). However, in the inhomogeneous model we must specify 
the behavior of ending spins in the large block limit. So we define the large block 
limit as when L — ^ cx), the two ending spins approach definite values, namely, 
Afoi S- and Ml,l+i — > 5'+. Then we realize from (j265p that as L — oo, 
fl^^j^ A(Z, A/j.j+i) — > Sifi. Therefore only the first term with Z = is left in (|275p . 
So that we need only to calculate the limiting Kq{{i): 

f / \ S— 



^l) (""l 



* \-./-+./ / t , },t 



(277) 
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Here both J_ and J+ take the hmiting values ^{S- — S+) and ^("S*- +5*+), respec- 
tively. 

Using transformation properties of the integrand in (|277p under SU{2), the 
KQ{Ct) integral is simplified and carried out as 

This expression states that {|G; J, f2)} is a set of eigenvectors of the density matrix 
as L ^ oo. Let us denote the density matrix in the limit by p^. Then (|278p leads 
to the result (see ([275ll ) 

pJG;^,f^) = ^^-^^l^|G;J,^7). (279) 
We find from (|279p that the limiting eigenvalue 

^-= (5_+i)W+i) ' "^^^ ^'''^ 

is independent of J. Any vector of the {S- + 1){S+ + l)-dimensional subspace 
spanned by the set {|G; J, f2)} is an eigenvector of p^^ with the same eigenvalue 
(S • Therefore p^ is proportional to a projector P(5'_+i)(5^+i): 

which is a generalization of (I230p to the inhomogeneous model. In the expression 
(|28ip S- and 5*+ denote the limiting spin values at the left and right boundary sites 
of the block, respectively. In addition, we also derive from the eigenvalues that the 
von Neumann entropy n coincides with the Renyi entropy Sb. and is equal to 
the saturated value In [(5_ + 1){S+ + 1)]. 



7. The One dimensional SU(n) Model 

In previous sections, we have discussed VBS states with spins in different repre- 
sentations of SU{2). Our discussion essentially exhausted all possible variations of 
SU (2) VBS states in 1-dimension with open boundary conditions. The AKLT model 
and the VBS ground state can be generalized to the SU (n) case as in I 28 | 29 | 68 |69| 
In this section we study the entanglement of the VBS state in the SU{n) version 
with open boundary conditions. Our treatment follows ^21. 



7.1. The Hamiltonian and the SU(n) VBS State 

Let us first define the model. Consider a 1-dimensional lattice with spins sitting on 
each site. What we mean by a 'spin' in our model is an adjoint representation of 
SU{n). Our spin chain consists of N adjoint representations of SU(n) in the bulk 
and fundamental and conjugate representations of SU{n) on the two boundaries. 
This construction corresponds to the spin-1 model in §|4] where we have the adjoint 
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representation of SU{2) (i.e. spin-1) in the bulk and the fundamental representation 
of SU{2) {i.e. spin-1/2) at two ends (the fundamental representation □ and its 
conjugate representation □ are equivalent for SU{2)). Unlike previous sections, in 
this section we shall reverse the order of discussion for better understanding. That 
is, we construct the SU{n) VBS state before writing down the Hamiltonian. The 
Hamiltonian will be constructed in such a way that the VBS state shall be the 
ground state. 

Now let us first construct an SU{n) VBS state which consists of N adjoint 
representations of SU{n) in the bulk and fundamental and conjugate representations 
of SU{n) on the boundary. First, we prepare sites k {k = 0, 1, ...,N) and k {k = 
1, 2, + 1) and arrange SU{n) singlets consisting of a fundamental (□) and its 
conjugate (□) representations as shown in FigureS] (see also^^^. We assign |j) (j = 



|0,0) |0,0/ |0,0) |0,0) |0,0) |0,0) |o,o) 




11 k k k+l k+l N N N+l 



Fig. 4. Construction of the SU (n) VBS state. White and black dots represent the SU (n) funda- 
mental and its conjugate representations, respectively. A solid line connecting two dots corresponds 
to the singlet state |0, 0) and a large circle denotes the projection onto the adjoint representation. 

0, 1, 2, n — 1) to the fundamental representation, while \j) {j ~ 0, 1, 2, n — 1) 
to the conjugate representation. \j) can be represented by the tensor product of 
{n- 1) \j)'s as 

J-—, E e^""--""l«2,...,a„) (282) 



where gJ"2,- -,a„ jg ^ totally antisymmetric tensor of rank n. Using \j) and \j), an 
SU{n) singlet state |0, 0) can be represented as a maximally entangled state: 



The above relation can be easily confirmed by inserting the resolution of the identity 
I = Yl'jZo and substituting (|282l) . Next, we prepare the adjoint representation 

of SU{n) by projecting the tensor product □ (g) □ onto an {n? — l)-dimensional 
subspace (the dimension of the adjoint representation of SU{n) is equal to the 
number of generators). This procedure corresponds to large circles in Figure ID In 
Figure [5ja) we visualize the decomposition rule □ ® □ = (singlet) © (adjoint). Then 
we have obtained the SU{n) adjoint representation at each composite site (k,k). 
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Henceforth we shall call this composite site k. Finally, we can represent the SU{n) 
generalized VBS state as 

N \ 



|VBS) = 




|0,0)oi|0,0)i2...|0,0) 



(284) 



where P^.^ is a projection operator onto an adjoint representation of SU{n). 

After construction of the VBS state ()284|) , we now write down the Hamiltonian 
along the same line as the SU (2) AKLT model: 

JV-l 

H=Y^ H{k, k + l) + H{0, 1) + H{N, N + 1), 



k=l 



H{k, k + l) = Y^ Cy(fc, k + 1) TTvik, k + 1), 



(285) 



where y is a Young tableau which is neither [n,n] nor [n, n — 1, 1]. Here we have 
assigned [ki, to the Young tableau Y, where kj is the number of boxes in the 

j*'^ column and Ai is the number of boxes in the first row. TTyik, fc + 1) is a projection 
operator which projects (adjoint) (adjoint) onto a representation characterized by 
Y and the coefficient Cyik, fc + 1) can be an arbitrary positive number. The reason 
why [n, n] and [n, n — 1, 1] are excluded from the sum is the following: Since □ at 
site k and □ at site fc + 1 have already formed a singlet in the ground state (12841) , the 
possible representations obtained from the decomposition of (adjoint) ® (adjoint) 
are restricted to [n, n] and [n, n— 1, 1] (as graphically shown in Figure[5jb)). i?(0, 1) 



□ 

(fundamental) 



> n-\ 



(conjugate) 



(b) 



V n 



(adjoint) 



(singlet) 



[«,«] [w,H— 1,1] 



Fig. 5. (a)The decomposition rule for the tensor product of □ ® (b)The Young tableaux 
corresponding to [n, n] and [n, n — 1, 1]. 



and H{N, N + 1) arc boundary terms which assure the uniqueness of the ground 
state of this Hamiltonian. i/(0, 1) and H{N, + 1) can be written in terms of 
the projection operators acting on the tensor products (fundamental) (8) (adjoint) 
and (conjugate) ig) (adjoint), respectively. By construction, the SU{n) VBS state 
(|284p is a zero-energy ground state of this Hamiltonian (|285p . We note here that 
another construction of the Hamiltonian by M. Greiter and S. Rachel I^Slig similar 
but slightly different from ours. 
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7.2. The Block Density Matrix 

Next, we consider the density matrix of a block subsystem of the VBS ground state 
|VBS) (|284p . To calculate the density matrix, it is more convenient to recast the 
chain of singlets |0, 0)oi|0, 0)i2 ■ ■ ■ |0, 0)^^^^^:^ in (|284p in a different form. Let us first 
look at a chain of two singlets |0, 0)oi|0, 0)i2- We can rewrite this product state as 

^ n — 1 n—1 

|0,0)oi|0,0)i2 = -^^ |Z,m)o2l^,-"i)li, (286) 

1=0 m=0 

where \m,n) is a basis of the maximally entangled state defined by 

\l,m) ^iUi^^®I)\0,0). (287) 

Here / is an n-dimensional identity matrix and Ui^m — Z'" (m, n = 0, 1, n—1) 
are generalized Pauli matrices, where the unitary operators X and Z act on |j) 
as X\j) = \j + l(mod.ri)) and Z\j) = with w = e^'^'/", respectively. One 

can easily show the relation (I286| by using the fact that |0, 0) is invariant under 
the action of (t//,™ ® Ui^-m), see^^. This procedure can be regarded as a multi- 
dimensional generalization of entanglement swapping. (The maximally entangled 
state basis ([M)) and a similar relation ([Ml have been used in ? I4.1I for the SU{2) 
spin-1 VBS state.) By repeatedly using the relation (|286p . we can generalize in a 
straightforward way to a chain of singlet states: 

|G,0)oi|0,0)i2...|0,0)^^ 

X! '■■ X! \h,-mi)ij^- ■ ■\lN,-mN)NN 

{h,mi) (/jv,"in) 

•(C/;i,mi • • • Ui^,m^ ^)|0, 0)ow+i^ (288) 

where {mk,nk) [k — 1,2, ...,N) runs from (0,0) to {n — l,n — 1). To obtain the 
ground state |VBS) from (|288p . we still have to make a projection onto the subspace 
of adjoint representation at each site k. Since the decomposition rule □ (g) □ = 
(singlet) ® (adjoint) and the fact that |0,0) is an SU{n) singlet, the vector space 
of the adjoint representation is spanned by |/, — m) ((/, to) ^ (0, 0)). Then the only 
thing to do is to omit the summation over {lk,mk) = (0,0) in (|288|) . As a result, 
the SU{n) generalized VBS state can be re-written as: 

\'^^^) = -(y^iym H \h,'rni)i^---\lN,-mN)nN 

7^(0,0) #(0,0) 

•(f/ii,mi • • • C^iiv,m« ^ 0)oN+-v (289) 

where we have already normalized |VBS) by the factor l/(n^ — 1)^/^. 

Now we calculate the density matrix of a block of contiguous spins of length L. 
Assume that the block starts from site k and stretches up to fc + L — 1, where k > 1 
and k + L ~ 1 < N. The reduced density matrix is obtained by taking the trace 
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over the sites j = 0, 1, k — 1 and j — k + L, N,N + I outside the block as 
PL = tri,...,fe-i,fe+L,...,Ar,o,lv+T [ I VBS) (VBS| ] 

7^(0,0) 5^(0,0) #(0,0) 

- E EE ••• E E 

*'i'/rn'(' 'iVn"n^ ''i^'"*' ('l, + fc-l;"'-L + fc-l) ('i, + fc_i.™i, + fc_i) 
5^(0,0) #(0,0) ^(0,0) #(0,0) ^(0,0) 

•|^fc, -TOfe)fefc(/fe,-774| • • • |/L+fc-l, -™L+fe-l)L + fe-lL+A;-l('L+fc-li-"^L+fc-ll 

■tra.W+i[{UiVU2<E>I)\0,0)„jr+i{0^(^\iUiV'U2<»T)^ , (290) 



where J7i — Ui^^mi ■ ■ ■ Ui^^umk^n U2 - UlL + k.rnL + k ■ • • Ul^^rnN^ 

y = Ui^,mk---Ui^^^_^,mr^+k_^ and V = Ui'^.m'^---Ui'^^^_^,„,'^^^_^. To rewrite the 
comphcated expression (|290l) . we use the following property of |0, 0): 

{S(g)T)\0,0) = {ST'(g>I)\0,0) = (/(8)rS'*)|0,0), (291) 

where S and T are n-dimensional unitary operations acting on \j) and |j), respec- 
tively, and the superscript t denotes the transposition. Using this property and the 
cyclic property of the trace, we can simplify the last part of ()290p as 

tro,lv+T [(U1VU2 (E> /)|0, 0)0^(0, 0\{UiV'U2 ® 1^] 
= K.-N+i [(Ui ® I){V (g, I){U2 ® /)|0,0)o^(0,0|(C/2 ® I)\V' (g> I)HUi (g> /)t] 
= tro.lv+T [(^ Ul}\0, 0)o^(0, 0|(/ ® U*^ {V ® /)t] 

= t%lv+T [(I ® U'2){V (g> /)|0, o)oiv+t(o> ® ® C^2)1 
= tro,iv+T[(^®^)|0,0)oiv+T(0,0|(V^'®/)t] . (292) 

Since (j292p does not depend on (/i, toi), • • • , {Ik-i, rrik-i) and {Ik-t-L^mk+L), • • • , 
(^AT, fcjv), we can rewrite the density matrix (j290p as 

Pl = ^ i)L EE"- E E h,-mk}-kk{l'k,-m'k\ 

/(0,0) #(0,0) #(0,0) ^(0,0) 

• • ■ Ih+k-l, -™-L+fc-l)r+irTL+fe_i(^L+fc-l> -ITT'L + k-ll 

•tro,lv+T[(^®^)|0,0)oiv+T(0,0|(F'®/)t] . (293) 

From the form of the density matrix in ([293|), we immediately notice that the density 
matrix does not depend on both the starting site k and the total length of the chain 
N. The same property for SU{2) VBS states has been discussed in S I4.2I . S j5^ and 
§ 16.21 (this property was proved in^ISlfQ]- gpin-l, in for Spin-S", and in for the 
inhomogeneous SU{2) models, respectively). We can regard the above result as an 
SU (n) generalization of those results. 
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7.3. Spectrum of the Density Matrix and Entropies 

Since the density matrix ()293|) is independent of both k and N, we can set 
N = L without loss of generahty. We can further reduce the original problem 
to that of the reduced density matrix of two ending spins (□ and □) using the 
Schmidt decomposition of a bipartite pure state (see the introduction §[1] for a brief 
discussion). Suppose that \B U E) is a bipartite pure state of a total system BU E. 
Then there exists orthonormal states \Bj) for the subsystem B, and orthonormal 
states \Ej) for E such that 

\BUE)=Y,V>^\Bj)®\E,}, (294) 
j 

where Xj{> 0) satisfy J2j — 1- The proof of the above theorem using the 
singular value decomposition can be found in §2.1.10 and §2.5 of From 
(|294p . one can immediately notice that the set of non- vanishing eigenvalues of 
Pb = tTE [\BU E)(BU E\] coincides with that of p^ = trs [\B U E)(B U E\]. 

Now we can reduce the eigenvalue-problem of pj^ to that of the density matrix 
for the two ending spins p^. This density matrix takes the following form: 

Pi = („2 ! E E iU^I)\0,0},—,{OMU^I)\ (295) 

(ii.™i) 

#(0,0) #(0,0) 

where U — C^ii.mi • • • Ui^^mL - To evaluate the eigenvalues of p^, it is convenient to 
formulate the action of (C/;,m /) as a transfer matrix. Let us first see the action 
of (f/i',m' (X" /) on a state \l,m): 

iJJi,,^,®I)\l,m) = (X''Z'"'X'Z™®/)|0,0) =cj"''|/ + r,7n + TO'), (296) 

where we have used the relation ZX = loXZ with uj = c^^'/". Using the above 
relation, we can prove that 

{Uv ^m' ® I)\l,m){l,ra\{Uv ^ra' ® I)^ = \l + l',m + m'){l + l',m + m'\. (297) 

Next, we assign the vector (0, 0, 1((/, to)*'' entry), 0, 0)* to the state \l,m){l,'m\. 
This one to one correspondence plays an essential role in our analysis. From this 
bijection, the operation 

T?x')#(o o)(^'',™' m\{Uv ^m' ® I)"^ can be written in terms of (n^ x v?)- 

dimensional matrix as 



i — 




/O 1 1 •• 




1 1 •• 


1 


1 1 •• 


1 


Viii.. 


0/ 
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This transfer matrix can be diagonalized by the foUowing unitary matrix: 



n 



(11 1 . 


1 


\ 


1 c • 




-1 


1 • 




-1) 


yl C"'"^ (^2(n"-l) . 







(299) 



where C = exp(27ri/n^). Then we can obtain the exphcit form of the reduced density 
matrix as 



PL 



1 



(n2 - l)L 
1 



T^(1,0,...,0)* 

C/,[diag(n2 - 1, -1, -l)]^t/t(l,0, 0)* 



(n2 - l)L 

^(l + (n2-lK(L))|0,0)„^(0,0| 



1 



^ (l-p„(L))|?,, 



(300) 



(i,m)7^(0,0) 



where we have used the relation, 1 + C*^ + C^'" + ' ' ' + C^"^ ^^'^ = (1 < fc < rP' — 1) 
and Pn[L) = { ^\ )^ . Substituting n = 2 into pOOl) . one can reproduce the result 

of the SU{2) Spin-1 VBS state obtained in §|1 (see also^S). 

Let us now start the evaluation of the von Neumann and the Renyi entropies of 
a block of L contiguous spins. First, we shall examine the von Neumann entropy of 
the block. From the Schmidt decomposition and the definition of the von Neumann 
entropy n [Pl] = S'v n [pi\ = -tri,2,...,L(PL InPi)' obtain 



Inn 



1 + (n2 - l)p„(L) 



ln(l + (n2 - l)p„(L)) 



-{n^-l)^-^4^Hl-Pn{L)) 



(301) 



with pn{L) ^ (;^)^. Similarly to the SU{2) VBS states HSHHESEZl and the XY 
spin chains in the gapped regime B4 | 45 | 22 | 2 ^^ n[Pl] is bounded by 2 Inn in the 
limit of large block sizes L — ?> cxd and approaches to this va lue exponentially fast in 
L. This is a partial proof of the conjecture proposed in^^, that the von Neumann 
entropy of a large block of spins in gapped spin chains shows saturation. Next we 
shall examine the Renyi entropy of our system. From the definition of the Renyi 
entropy Sjx{a) = ^ In [tr(p2)] (a 7^ 1 and a > 0), 



1-a 



ln(Ao,o(i)" + (ri'-l)Ai,™^o,o(i)"), 



(302) 
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where 

( ^il + in^-l)PniL)), (;,m) = (0,0) 
\m{L) = <^ (303) 

I ^il-Pn{L)) (1,171)^(0,0). 

7.4. The Density Matrix as a Projector 

We have obtained the spectrum and derived the entropies of the density matrix in 
S I7.3I . Our treatment avoided exphcit construction of the eigenvectors by reducing 
the problem to the density matrix of the two ending spins. In this section, we shah 
show that Theorem [2] in i^ l3.3l is also valid for our SU(n) version VBS state, i.e. 
The eigenvectors of with non-zero eigenvalues are degenerate ground states of 
the block Hamiltonian (The block Hamiltonian is defined similarly as in the SU(2) 
model, see p04[) below). Indeed, the proof of Theorem [5] is completely applicable 
in our SU (n) case because it only relies on the definition of the density matrix and 
the fact that the VBS state has no projection on any of the subspaces specified in 
the Hamiltonian. 

The SU(n) block Hamiltonian, defined along the same line as in the SU(2) 
model, is the sum of interaction terms within the block: 

L-l 

Hb ^J2H(k,k + l). (304) 
fe=i 

The number of degrees of freedom in the block subsystem can be counted from the 
von Neumann entropy as . This number coincides with the number of degenerate 
ground states of the block Hamiltonian Hb- i.e. 

D = deg = n^, (305) 

where D is the number of non-zero eigenvalues (dimension of the support) of the 
density matrix and deg is the number of degenerate ground states of the block 
Hamiltonian. A basis of the degenerate ground states can be constructed as follows: 

\YBS; p, q) = Cp^q ^ ••• ^ -TOi)n • • • -mL_i)— 

#(0,0) /(0,0) 

•Pll {(Up.^gUi,,^,...Ui,_,.„,,-, ® I)\0,0)l-l) , (306) 

where Cp^q is a normalization factor and Pj^i is a projector onto the adjoint rep- 
resentation of SU(n). This set of states (I306P can be called the degenerate VBS 
states. Any linear combination of (|306p is apparently a ground state of Hb (|304p . 
The graphical representation of the construction of this state is shown in Figured 
The following orthogonality relation of these generate VBS states (I306P holds: 



(VBS;p,g|VBS;r,s) = Clq(n^ - l)^,5p,,,5,,,A_p,_,(L) 



(307) 
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\VBS;p,q) 




1 1 



L L 



Fig. 6. Construction of the degenerate VBS states |VBS;p, 5) for the block. A white (black) dot 
represents the SU (n) fundamental (conjugate) representation. A large circle denotes the projection 
onto the adjoint representation. The dashed line corresponds to the state |VBS;p, g). 



where the subscripts — p and —q are modulo n. This can be shown as follows: 
(VBS;p,q|VBS;r,s) 



— Cp^qCr^s • • • 



7^(0,0) 5^(0,0) 

LL i,mL-i)^PLL(tA-,sC^;i,mi---f^ii,_i,mi_i ® ^) LL 

#(0,0) #(0,0) #(0,0) 
LL{0,0\iUl,^rm-Ul^_,^m^_,® lYiUp^q^ I)^l,m)LL 
■LL{l,m\{Ur,s ® /)(f//i,mi-f//i_i,m^_i I)\0,0)lL 
— Cp,qCr^s{n-'^ ^ 1)"^ 

• Xi,,,n'iL~l)LL{p + l',q + m'\il-\0,0)LL{0,0\)\r + l\s-^m')LL 



{l',m') 



(308) 



Here we have recalled (|295p and used the relation 



#(0,0) #(0,0) 



((7,,,,„,...t/z^_,,™^_, ®/)|0,0)i£(0,0|((7,,,„,...[/z^_,,™^_, 
= (n^ - 1)^-Vrri = (^i' - 1)^"' E ^'.™(^ - l)|?,™)L£(/,m|. (309) 

(i,m) 

The explicit form of the normalization factors Cp^q are given by Cj 



-'PA 



1 / \J {rfi — l)-^A_p__g(i). Now, we could write in terms of this basis of degenerate 
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VBS states. By the original definition, 
PL =tro— [|VBS)(VBS|] 



E 


E 






#(0,0) 


#(0,0) 


mL-i){l'L_ 





1^2 „ 1)L ^ ^ ^ 

#(0,0) 5^(0,0) 

|/i,-mi)(/^,-m'i| • • • 

•PLL(t/p,«C^ii,mi---C^i.-i,m.-i®^)|0,0)i£(0,0|(C/p^,t/r^,„;...C/r^_^,W,_,®/)^PLL 

(310) 

Then by comparing with p06l) . we obtain 

Pl^Y. A-p.-g(i)|VBS;p,g)(VBS;p,<?|, (311) 

(P,9) 

where \^p-q{— \n-p,n~q) was defined in p03p . Therefore we conclude that the 
density matrix of a block of L contiguous spins in the SU{n) VBS state is 
completely characterized by the degenerate ground states { |VBS;p, q) } of the block 
Hamiltonian Hb ([304ll . 

7.5. T/ie Large Block Limit 

Now we consider the large block limit, i.e. L — cx). In this case, Pn{L) — and Aj.m 
(|303p become degenerate. So that great simplification occurs in the expressions of 
entropies as L — > oo: 

S'v N = In 71^, 
SK{a) = -^In f fijV + (n^ - 1) ("^ 

^ ^_ln(„2)i-" 

1 — a 

= ln7i2. (312) 

We notice that the Renyi entropy is independent of a and coincides with the von 
Neumann entropy in the large block limit. The same saturation behavior was ob- 
served in all our SU{2) cases in § 14.51 . 5 15. 81 . S I6.6I . This means that the density 
matrix of a large block is proportional to a n^-dimensional identity matrix. In other 
words, a sufficiently large block of neighboring spins in our SU (n) VBS ground state 
is maximally entangled with the rest of the chain. 

In the limit of large block sizes, i.e. L — > oo, pj^ can be written as 

Pl^Poo = AElVBS;p,g)(VBS;p,g| = ^P„2, L ^ ^. (313) 

(P,9) 

The limiting density matrix is proportional to a projector P„2 which projects on 
the n^-dimensional subspace spanned by the degenerate ground states of the block 
Hamiltonian. This structure of the limiting density matrix is also a generalization 
of the corresponding results for SU(2) density matrices in §|43], ^ \5.8\ . ^ \6.6\ . 
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8. Summary 

We have studied the entanglement in the VBS ground state of the AKLT model in 
this review. The AKLT model is formulated on an arbitrary connected graph or a 
lattice. The Hamiltonian ([20|) . (|28l) is a sum of projectors which describe interactions 
between nearest neighbors. The condition of uniqueness of the ground state relates 
the spin value at each vertex (site) with multiplicity numbers associated with edges 
incident to the vertex (bonds connected to the site), see (HJ), (1511) . The unique 
ground state is known as the Valence-Bond-Solid state (|26|). ((38|). 

To study the entanglement, the graph (lattice) is divided into two parts: the block 
and the environment. We investigate the density matrix pg of the block and show 
that it has many zero eigenvalues. We describe the subspace (called the support) 
of eigenvectors of with non-zero eigenvalues. It has been proved (see Theorem 
[2] in § 13.31 ) that this subspace is the degenerate ground space of some Hamiltonian 
which is called the block Hamiltonian Hb (|68p . The block Hamiltonian is a part of 
the original AKLT Hamiltonian describing interactions of spins inside of the block. 

The entanglement can be measured by the von Neumann entropy or the Renyi 
entropy of the density matrix pg. Most eigenvalues of pg vanish and have no 
contribution to the entanglement entropies. The density matrix takes the form of a 
projector on the ground space of Hb multiplied by another matrix (see also . 




Fig. 7. A 2-dimensional Cayley tree. Each dot represents a spin-35/2 in the bulk and a spin-5/2 
on the boundary. Each solid line represents the bond connecting a pair of interacting spins. This 
tree structure has no loop. 

A complete analysis of the density matrix for a variety of 1-dimensional AKLT 
models (including the SU{n) generalization) has been presented. The block density 
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matrix pj^ for a subsystem of L contiKuo us bulk spins h as been diagonalized with 



the density matrix pg is changed to pj^ for these 1-dimensional models as empha- 
sizing the dependence on the size L of the block.) We find that in all these cases 
the support coincides with the ground space, so their dimensions are equal D = deg. 
In the large block limit L all non-zero eigenvalues become the same and 

the density matrix is proportional to a projector (|129p . (I230p . (I28ip . (|313p . The 
von Neumann entropy equals the Renyi entropy and both take the saturated value 
S'v N = '5'r = In = ln(c?eg). 

Moreover, it turns out that the block Hamiltonian Hb defines the density matrix 
Pj^ completely in the large block limit L — ^ cx). The zero-energy ground states of 
the block Hamiltonian Hb span the subspace that the density matrix pj^ projects 
onto. So that can be represented as the zero-temperature limit of the canonical 
ensemble density matrix defined by Hb- 



In the zero-temperature limit, contributions from excited states of Hb all vanish 
and the right hand side of pi4p turns into a projector onto the ground states of 
the block Hamiltonian. 

For more complicated graphs or lattices, non-zero eigenvalues of the density 
matrix are still unknown. One open problem is to calculate those eigenvalues. One 
may start with the Cayley tree (also known as the Bethe tree), see Figure [71 The 
picture shows a Cayley tree with each bulk vertex connected to three edges. The 
uniqueness condition (|30|) requires that we shall have spin-3S'/2 in the bulk and 
spin-S'/2 on the boundary. A symmetric block subsystem consists of spins enclosed 
by a circle centered at the center of the tree. The degeneracy of ground states of the 
block Hamiltonian is (^-l- 1)^^^, where Nqb is the number of sites on the boundary 
of the block. An exact explicit expression for the non-zero eigenvalues (finite block) 
is expected because there is no loop. It i s also interesting to study the large block 
limit. In all known examples ^ISESEdIZZI -^^here the density matrix of a large block 
has been calculated, all non-zero eigenvalues approach the same value 1/D = l/deg. 
So that the entanglement entropies are saturated, i.e. S'v n = ■S'r = ln£) = In (deg). 
Therefore the density matrix of a large block is proportional to a projector on the 
ground space of the block Hamiltonian, i.e. Pb — J)Pd — '^^deg- However, this 
might not be the case for the Cayley tree. According to the area law, we expect 
that in the large block limit (a circular block centered at the center of the tree), the 
entropy be proportional to the length of the boundary with some coefficient a, i.e. 
S'v N = oiNqb- It is interesting to calculate a and we expect that it will be smaller 
than ln(5 -I- 1). Another open problem is the generalization of the VBS state to 
other Lie algebras beyond SU{n) and the study of the entanglement. 



non-zero eigenvalues calculated 





L — > oo. 



(314) 



"As L — > oo, the size of the whole spin chain also goes to infinity. 
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It is also important to calculate non-zero eigenvalues of pg for graphs with loops. 
For example, consider a basic AKLT model defined on the 2-dimensional hexagonal 




Fig. 8. The basic model on a 2-dimensional hexagonal lattice. Each spin-3/2 in the bulk is rep- 
resented by three small dots representing spin-l/2's enclosed by a solid circle (a symmetrization). 
The solid lines antisymmetrize the connected two spin-1/2 states. The block is given by a dashed 
large circle. Each line cut by the dashed circle results in a free spin-1/2 on the boundary of the 
block. 

lattice, see Figure [D The basic model has spin-3/2 in the bulk and the block is a 
large circle. Note that each spin on the boundary of the block contributes effectively 
a free spin-1/2. Then according to Kastura's formula ((75)) . the degeneracy of the 
ground states of the block Hamiltonian is 2^^^ , where Nqb denotes the number of 
site (spins) on the boundary of the block. There is at most 2^^^ number of non-zero 
eigenvalues of the density matrix. The entropy takes a saturated value if all these 
eigenvalues are equal in the large block limit. Similar to the case of the Cayley tree, 
in the large block limit the entropy should be proportional to the size (area) of the 
boundary, i.e. S'v n = ctNgg. We expect that the coefficient a be greater than 
but smaller than In 2. The value Nqb In 2 is an upper bound of the von Neumann 
entropy because this is the logarithm of the dimension of the Hilbert space (number 
of non-zero eigenvalues of the density matrix) . 
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